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RESUMO 


Neste artigo descrevem-se os 
mecanismos de aceleração 
de partículas carregadas em 
plasmas criados por lasers ul- 
tra intensos e apresenta-se O 
princípio físico subjacente a 
tecnologia dos sistemas laser 
TeraWatt, Descreve-se ainda o 
sistema TeraWatt a instalar no 
IST para a realização de expe- 
riências de aceleração de 
fotões, geração de ondas 
plasma electrónicas de grande 
amplitude e aceleração de 
electrões. 


Abstract 


In this paper, we describe the 
mechanisms of plasma particle 
acceleration in laser-plasma in- 
teraction schemes and we pre- 
sent the technology of Tera- 
Watt laser system, and its 
underlying physical principle. 
The TeraWatt system to be 
installed at IST is described, 
This system will be the basis 
for photon acceleration experi- 
ments, laser excitation of 
large-amplitude plasma waves, 
and electron acceleration ex- 
periments to be undertaken at 
IST. 


1. INTRODUÇÃO 


Nos últimos sessenta anos, os 
avanços na tecnologia dos acelera- 
dores de partículas têm sido consi- 
deráveis. Acelerando partículas car- 
regadas até velocidades próximas 
da velocidade da luz e fazendo-as 
colidir, tem permitido avanços signi- 
ficativos no conhecimento das for- 
ças fundamentais da natureza, da 
origem do universo e da estrutura 
da matéria. O aumento da energia 
das partículas aceleradas tem per- 
mitido estudar as reacções nuclea- 
res no interior das estrelas gigan- 
tes, o núcleo das estrelas de 
neutrões e, com as máquinas mais 
recentes, as condições do universo 
até ao instante t-=10º s após a sua 
criação tal como o descreve o mo- 
delo do big-bang. No entanto, os 
aceleradores actualmente em cons- 
trução ou em projecto já começam 
a aproximar-se dos limites tecnoló- 
gicos actuais e os custos orçamen- 
tais atingem valores incompatíveis 
com o estado actual da economia 
mundial. De facto, o projecto de 
construção do maior acelerador de 
partículas do mundo (SSC — Super- 
conducting Supercollider) foi aban- 
donado devido aos elevados custos 
envolvidos na sua construção [1]. 
O SSC estaria apetrechado com 
um anel de aceleração de 87 km 
que possibitaria a aceleração de 
partículas até aos 40 TeV. O ele- 
vado investimento necessário é, em 
parte, devido ao próprio conceito de 
acelerador que se pretendia imple- 
mentar: as partículas são guiadas 
por campos magnéticos ao longo 
do túnel de aceleração e acelera- 
das por campos eléctricos intensos. 
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Para reduzir as dimensões de um 
tal acelerador seria necessário au- 
mentar consideravelmente a intensi- 
dade dos campos magnéticos e 
eléctricos necessários para guiar e 
acelerar as partículas carregadas. 
Contudo, os limites tecnológicos 
não permitem encarar esta hipótese 
com esperança. As forças magnéti- 
cas nos aceleradores de partículas 
convencionais já são da ordem de 
grandeza das forças que preservam 
a estrutura do material magnético. 
Por outro lado, os campos eléctri- 
cos utilizados atingem valores pró- 
ximos dos campos de ligação dos 
electrões nos átomos, o que signi- 
fica que campos eléctricos mais 
elevados arrancariam os electrões 
dos átomos que constituem a estru- 
tura dos aceleradores de partículas 
destruindo esses materiais. Desta 
forma, tal como anteriormente se 
passou dos aceleradores electrostá- 
ticos para os betatrões e ciclotrões, 
e destes até às máquinas actuais 
(os aneis de colisão), é necessário 
agora explorar novos conceitos de 
aceleradores baseados em proces- 
sos físicos novos, como é o caso 
dos aceleradores de partículas a 
plasma. 


2. ACELERADORES A PLASMA 


O plasma é um estado da matéria 
onde os electrões se encontram 
dissociados dos átomos que consti- 
tuem o meio. Apesar deste facto, o 
meio é globalmente neutro sendo 
possível criar zonas deslocalizadas 
de carga d correspondendo a zonas 
de carga negativa positiva i.e. zo- 
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nas em que a densidade electró- 
nica é superior ou inferior à densi- 
dade iónica. Nos aceleradores de 
particulas a plasma, os campos 
eléctricos presentes não estão limi- 
tados pelas energias de ligação dos 
electrões aos átomos. Estes acele- 
radores baseiam-se no princípio de 
que as partículas carregadas são 
aceleradas por campos eléctricos 
intensos, no interior do plasma. 
Dado que o plasma corresponde a 
um estado da matéria em que as 
energias dos electrões são suficien- 
tes para estes estarem completa- 
mente dissociados dos átomos a 
que pertenciam, e portanto o plas- 
ma já se encontra ionizado, os 
campos eléctricos que o plasma 
pode sustentar são bastante mais 
elevados já não existe a possibli- 
dade de dissociação electrónica. 
Esta é a vantagem decisiva deste 
novo conceito de acelerador. Com 
efeito, teoricamente, um plasma 
pode sustentar campos eléctricos 
10000 vezes mais fortes (10º-10º 
V/cm) do que aqueles presentes 
nos aceleradores convencionais, 
para valores típicos da densidade 
electrónica. Isto significa que se 
campos eléctricos desta ordem pu- 
derem ser produzidos num plasma 
ao longo de uma distância conside- 
ravel, um acelerador a plasma com 
algumas centenas de metros de 
comprimento será capaz de produ- 
zir resultados semelhantes as 
previstas para o SSC. Vejamos re- 
sumidamente os antecedentes his- 
tóricos dos actuais conceitos de 
aceleradores a plasma. A acelera- 
ção de partículas carregadas por 
campos magnéticos movendo-se no 
espaço e por ondas electromagnéti- 
cas foi inicialmente considerada por 
Fermi [2] e McMillan [3] quando nos 
anos 40 e 50 se ocuparam da aná- 
lise dos raios cósmicos. O mesmo 
problema foi posteriormente tratado 
por Sagdeev [4], considerando 
agora o plasma e não o vazio como 
o meio no qual se movem as partí- 
culas e se geram os campos. Esta 
inclusão do plasma na problemática 
da aceleração foi particularmente 
importante, dado que o plasma per- 
mite a excitação de campos eléctri- 
cos extremamente elevados. Com 


efeito, uma pequena perturbação é 
suficiente para estabelecer no 
plasma uma oscilação da densi- 
dade electrónica, de frequência «o, 
vn, em que w, se designa por fre- 
quência plasma e ne é a densidade 
electrônica. Estas oscilações ou on- 
das do plasma apresentam uma 
propriedade fundamental: a sua fre- 
quência é praticamente indepen- 
dente do comprimento de onda. Da- 
qui resulta que a velocidade de 
fase das oscilações do plasma 
pode tomar qualquer valores arbi- 
trariamente elevados desde que Oo 
comprimento de onda seja esco- 
lhido de forma apropriada. Por ou- 
tro lado, a velocidade de grupo das 
oscilações plasma-electrónicas é 
nula e desta forma, assim que a 
onda plasma electrônica é estabele- 
cida no plasma, a sua amplitude 
mantém-se estacionária até ser 
amortecida ou destruída por algum 
processo não linear. Observe-se 
ainda que o conceito de aceleração 
de partículas que está subjacente à 
física das altas energias, afasta-se 
bastante do conceito intuitivo de va- 
riação da velocidade da partícula 
ao longo do tempo. Nos acelerado- 
res, as partículas movem-se a velo- 
cidades próximas da velocidade da 
luz e quando absorvem energia por 
interacção com os campos, são 
aceleradas na medida em que a 
sua massa efectiva aumenta, de 
acordo com a teoria da relatividade 
restrita. A velocidade das partículas 
aceleradas apemas varia ligeira- 
mente uma vez que está perto do 
limite físico que é a velocidade da 
luz, Desta forma quando se pre- 
tende acelerar particulas relativistas 
é necessário que as oscilações 
plasma electrônicas se movam com 
velocidades próximas de c de modo 
a que os electrões a acelerar acom- 
panhem as ondas plasma electróni- 
cas. 


3. BEAT-WAVE E WAKEFIELD 


Apesar dos plasmas permitirem a 
existência de ondas com velocida- 
des de fase arbitrária, e com cam- 
pos eléctricos intensos, resta resol- 
ver o problema da excitação dessas 


ondas. Dois métodos foram propos- 
tos. O primeiro método, denomi- 
nado por beat-wave, basela-se na 
co-propagação através do plasma 
de dois feixes laser de frequências 
w, é w, À diferença de frequência 
D,-), é escolhida de modo a coinci- 
dir com a frequência plasma: w,- 
-D=0Wp. À velocidade de fase das 
ondas plasma electrónicas excita- 
das por este processo é igual à ve- 
locidade de grupo do batimento dos 
dois feixes laser que, para frequên- 
cias w, e w, muito superiores à 
quando comparadas com a fre- 
quência plasma q, é muito próxima 
de c. Este método foi proposto por 
Tajima e Dawson em 1979 [5]. 

O segundo método baseia-se na 
propagação de um impulso laser 
de duração n= EE através de um 
plasma. Nesta situação, o laser 
deixa atrás de si uma espécie de 
rasto que é essencialmente uma 
oscilação plasma electrónica, de ve- 
locidade de fase igual à velocidade 
de grupo do impulso laser. A ampli- 
tude da onda plasma electrónica 
será proporcional à intensidade do 
impulso laser. Este método de- 
signa-se por aceleração de parti- 
culas por wake field [6]. 

As dificuldades experimentais asso- 
ciadas à aceleração de partículas 
por beat wave são óbvias. Em pri- 
meiro lugar, é necessário garantir 
que a diferença de frequência entre 
os dois feixes laser seja igual à fre- 
quência plasma electrónica do 
plasma onde os dois feixes se pro- 
pagam. Por outro lado, a duração 
dos impulsos laser tem que ser su- 
ficientemente pequena de forma a 
não perturbar o movimento dos 
iões do plasma. O crescimento das 
ondas plasma electrónicas por beat 
wave é um processo não linear que 
é ressonante quando w,-0,=0, [7]. 
Quando os dois feixes laser se 
propagam no plasma, actuam so- 
bre os electrões forças não lineares 
jixB; e j;xB, que dependem da 
frequência plasma. Estas forças 
são paralelas e anti-paralelas à di- 
recção de propagação dos impulsos 


laser e são responsáveis pelo cres- 
cimento ressonante das ondas 
plasma electrónicas com o número 
de onda k=k,-k, e frequência 0,. 
Este crescimento será tanto mais 
eficiente quanto menor for a dife- 
rença entre a frequência plasma 
electrónica do meio e a diferença 
de frequência w,-w, Apesar de ser 
possível criar plasmas com densi- 
dades bem definidas, uma pequena 
diferença é suficiente para que não 
se atinja a situação de ressonância 
e para que a geração da onda 
plasma electrónica de grande am- 
plitude não se observe. O outro 
problema experimental já referido 
prende-se com a duração dos im- 
pulsos laser. Com efeito, se a dura- 
ção dos impulsos laser for suficien- 
temente pequena (Tr < O.ins) o 
campo eléctrico do laser arrastará 
somente os electrões e os iões, de 
massa muito superior à dos elec- 
trões, sofrerão uma aceleração bas- 
tante mais pequena durante um in- 
tervalo de tempo reduzido; 
consequentemente, o movimento 
ióúnico será desprezável, Caso esta 
situação não se verifique, o movi- 
mento iónico é responsável pela ini- 
bição do crescimento da onda 
plasma electrónica através do cres- 
cimento de instabilidades iónicas 
do desenvolvimento de turbulência 
plasma. 


Apesar destes condicionantes, os 
resultados experimentais relativos a 
este método são bastante animado- 
res. Uma experiência conduzida por 
C.Joshi em U.C.L.A. [8] revelou a 
presença de campos eléctricos da 
ordem dos 10ºV/cm em ondas elec- 
trostáticas propagando-se com ve- 
locidades de fase próximas de c. 
Nesta experiência utilizou-se um la- 
ser de CO, de 70 J num impulso de 
300 ps que pode ser divido em dois 
feixes de comprimentos de onda 
10.29 um e 10.59 um de energias 
10 J e 60 J, respectivamente. 
O plasma é criado através de uma 
descarga de arco numa câmara 
contendo Hidrogénio à pressão de 
110-200 mTorr. Os dois feixes fo- 
ram opticamente misturados no 
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plasma criando ondas plasma de 
grande amplitude. As caracteristi- 
cas das ondas plasma electrónicas 
criadas foram analizadas através de 
diagnósticos ópticos (Difusão de 
Thomson). Injectaram-se electrões 
de energia 2.1 MeV que foram ace- 
lerados até energias próximas dos 
30 MeV, ao longo de um 1cm. 
O conceito de aceleração de parti- 
culas por beat-wave não é exclu- 
sivo da física dos aceleradores a 
plasmas; para regimes de energia 
mais baixos utiliza-se o mesmo 
conceito na geração de corrente 
não indutiva em plasmas de fusão 
nuclear [9]. A aceleração de parti- 
culas por wake field é, aparente- 
mente, mais fácil de implementar 
do ponto de vista experimental em- 
bora não tenha ainda sido demons- 
trado experimentalmente. O fenó- 
meno físico que explica a geração 
de ondas plasma electrónicas por 
wake field é extremamente simples: 
um feixe laser que se propaga num 
plasma empurra os electrões que 
encontra ao longo do seu percurso. 
Após atravessar essa região, 
forma-se uma zona de carga espa- 
cial positiva com a dimensão espa- 
cial do impulso laser e consequen- 
temente cria-se um campo eléctrico 
que arrasta os electrões para a 
zona de carga positiva e desenca- 
deia a oscilação plasma electrónica. 
Esta onda electrostática move-se 
com a velocidade de grupo do im- 
pulso laser e o seu crescimento 
não não é condicionado por ne- 
nhum fenómeno ressonante. Ape- 
nas é necessário garantir que o in- 
verso da duração do impulso laser 
seja da ordem da frequência plas- 
ma electrónica e que a energia do 
impulso laser seja elevada. Este 
facto permite maior flexibilidade 
quanto às características do plasma 
onde serão geradas estas oscila- 
ções, ao contrário do que sucede 
no beat-wave. No entanto, levan- 
tam-se problemas técnicos impor- 
tantes relativamente às característi- 
cas do impulso laser (duração 
inferior a 1 ps e energia da ordem 
das centenas de md) a utilizar nas 


experiências de wake field. As ex- 
periências com estes lasers inicia- 
ram-se recentemente e ainda fo- 
ram conclusivas: a presença do 
wake field ainda não foi detectada. 
Actualmente, preparam-se expe- 
riências de aceleração de particulas 
por plasmas em U.C.L.A. e La- 
wrence Livermore Laboratory 
(U.S.A). LULI-Ecole Polytechnique 
(França), Rutherford Appleton Labo- 
ratory (UK), Universidade de Osaka 
(Japão). Os trabalhos preparatórios 
para a experiência portuguesa de 
aceleração de particulas por wake 
field foram iniciados em Janeiro de 
1994 no Centro de Electrodinâmica 
do IST. 


4. O SURFATRÃO 


A aceleração de partículas por on- 
das plasma electrónicas de grande 
amplitude apresenta dois aspectos 
tecnicos fundamentais. O primeiro 
aspecto prende-se com o próprio 
crescimento das ondas plasma 
electrónicas e com a sua amplitude 
máxima. De facto, existe um limite 
maximo para a energia que uma 
onda plasma electrónica pode 
transmitir a um electrão. Esse limite 
é estabelecido pelo campo eléctrico 
máximo que pode ser produzido no 
plasma. Este limite corresponde às 
oscilações máximas de densidade 
electrónica que é determinada pe- 
las condições em que a oscilação é 
destruída por processos não linea- 
res ( wave breaking) [10]. Esta des- 
truição é semelhante à que ocorre 
nas ondas do mar quando a crista 
da onda desaba na base da onda 
— rebentação. Teoricamente, pre- 
vê-se que a amplitude máxima, 
para a qual não existe rebentação 
da onda plasma electrónica, é pro- 
porcional à raiz quadrada da densi- 
dade electrónica, o que nos permite 
obter valores do campo eléctrico da 
ordem das dezenas de GV/m. 

O segundo aspecto relaciona-se 
com a dinâmica do electrão no 
campo criado pela onda. Verifica-se 


intuitivamente, através de uma ana- 
logia para o caso clássico de uma 
partícula sujeita a um potencial si- 
nusoidal [11], que o tempo de inte- 
racção entre a partícula carregada 
e a onda electrostática é reduzido 
se a velocidade do electrão for su- 
perior à velocidade da onda plasma 
electrónica. Apesar destas velocida- 
des serem bastante próximas, 
o electrão ultrapassará a onda 
plasma e o tempo de interacção 
assim como a energia cedida pela 
onda ao electrão serão reduzidos. 
A forma de evitar esta situação 
consiste em prolongar o tempo de 
interacção do electrão com a onda. 
Tal como um surfista se move com 
um determinado ângulo em relação 
à direcção de propagação da onda 
em que se desloca de forma a au- 
mentar o tempo de interacção com 
a onda e a aumentar a sua veloci- 
dade, também um electrão mo- 
vendo-se obliquamente em relação 
à direcção de propagação da onda 
estará mais tempo sob a influência 
dessa onda e, portanto, a energia 
de interacção será mais elevada. 
Neste caso, a velocidade dos elec- 
trões tem que ser maior que a velo- 
cidade da perturbação de forma a 
que os electrões acompanhem a 
onda. Com o movimento oblíquo, o 
desfasamento entre o electrão e a 
onda, que aumentaria continua- 
mente, devido à diferença de velo- 
cidades, é eliminado. A forma de 
impór aos electrões um movimento 
obliquo em relação a direcção de 
propagação da onda consiste em 
aplicar um campo magnético verti- 
cal e constante ao sistema:o elec- 
trão fica sujeito a uma força de 
deriva E xB cuja direcção é per- 
pendicular à direcção de propaga- 
ção da onda electrostática (Figura 
1). Esta técnica, proposta por Kat- 
souleas e Dawson em 1983 [12], 
designa-se por Surfatrão, por ana- 
logia óbvia com a modalidade des- 
portiva que serviu de inspiração ao 
conceito. 

A situação mais geral de uma par- 
ticula carregada propagando-se no 
campo de uma onda electrostática 


sob acção de um campo magné- 
tico constante é bem conhecida 
[13] e verifica-se que, em determi- 
nadas condições, a aceleração de 
partículas passa a ser estocástica. 
Neste regime, um feixe de elec- 
trões injectado com uma determi- 
nada velocidade sofrerá um alarga- 
mento do seu espectro de energia. 
Embora exista aceleração das par- 
tículas carregadas nesta situação, 
a sua eficiência é reduzida e as 
energias alcançadas consideravel- 
mente mais baixas do que na ace- 
leração regular. 


5. LASERS DE TERAWATT 


A aceleração de partículas por 
plasma desenvolveu-se nos últimos 
anos de maneira acentuada devido 
aos avanços tecnológicos no domi- 
nio dos lasers de impulsos ultra 
curtos e intensidades elevadas. 
Desde 1985, é possível conceber 
lasers com durações inferiores ao 
picosegundo e energias da ordem 
do Joule [14]. A estes valores cor- 
respondem potências superiores ao 
TeraWatt. As energias disponíveis 
e a sua duração tornam estes siste- 
mas ideais para as experiências de 
aceleração de partículas em plas- 
mas e, de forma geral, para o es- 
tudo dos fenômenos de interacção 
de campos electromagnéticos ultra- 
"intensos com a matéria. Dada a in- 
tensidade dos campos eléctricos 
gerados pelo laser de TeraWatt, o 
estudo dos fenómenos físicos de in- 
teracção laser-matéria neste regime 
requer um tratamento relativista e 
não linear. Por outro lado, esses 
mesmos lasers permitem estudos 
em áreas tão distintas como a diná- 
mica das reacções químicas ultra- 
-rápidas, geração de harmónicas 
em gases, fontes de raios-x para a 
biologia e medicina, espectroscopia 
rápida, diagnósticos de fusão nu- 
clear, caos quântico, óptica quân- 
tica, geração de luz branca em li- 
quidos, etc. 


A produção de impulsos de duração 
inferior ao picosegundo está asso- 
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ciada a uma técnica proposta em 
1985 por Moorou e designada por 
CPA ( Chirped Pulse Amplification). 
O principio físico subjacente a esta 
técnica é simples; se pretendemos 
impulsos extremamente curtos en- 
tão a largura espectral desses im- 
pulsos deve ser bastante elevada 
(AmaAt < 1). Por outro lado, os am- 
plificadores não suportam potências 
excessivamente elevadas. A técnica 
CPA resulta, em traços gerais, des- 
tas limitações: o impulso inicial do 
laser é alargado antes de ser ampli- 
ficado de modo a não danificar os 
amplificadores. Seguidamente, o 
impulso é comprimido até durações 
próximas daquelas com que foi 
injectado na fase de alargamento. 
À saída do sistema, a potência 
atinje valores superiores ao Tera- 
Watt: com tecnologia Ti:S as dura- 
ções podem ir até às dezenas de 
femtosegundo e dezenas de mili- 
joule enquanto que com osciladores 
de Nd atinjem-se energias de 10J 
com impulsos de 1 picosegundo. 


O sistema CPA-TeraWatt de Lisboa 
(Figura 2), em fase de projecto e 
montagem, é um exemplo desta 
técnica. O impulso de um oscilador 
Nd:YLF bombeado por um laser de 
diodo é alargado por meio de redes 
de difracção e seguidamente ampli- 
ficado num amplificador regenera- 
tivo de Ti:S e num amplificador de 
Nd:Vidro. Após a fase de amplifica- 
ção é comprimido com a ajuda de 
redes de difracção, resultando num 
impulso laser de energia 1J e dura- 
ção inferior a 1ps, com um compri- 
mento de onda da ordem de 1 mi- 
cron (infra vermelho próximo). 


Este sistema constituirá o equipa- 
mento base para a experiência por- 
tuguesa de aceleração de parti- 
culas em plasmas que será 
desenvolvida em colaboração com 
LULI-École Polytechnique e RAL. 
Os objectivos do projecto português 
consistem em instalar este sistema 
laser em Portugal, observar 
experimentalmente a presença do 
wake field e desenvolver seguida- 
mente uma experiência de acelera- 


ção de electrões. Pretendem-se 
ainda, paralelamente, avançar com 
experiências de aceleração de fo- 
tões com o objectivo de desenvol- 
ver novas fontes de radiamento. Tal 
como na aceleração de partículas, 
a interacção de um fotão com uma 
perturbação da densidade electró- 
nica do plasma que se desloca ra- 
pidamente no espaço leva a altera- 
ção da energia dos fotões que 
interagem com essa perturbação. 
Esta alteração tem como conse- 
quência uma variação da frequên- 
cia do fotão dependente da veloci- 
dade da perturbação [15]. Desta 
forma é possível conceber fontes 
de radiamento sintonizáveis, desde 
as micro-ondas até aos raios-X, As 
potencialidades do sistema laser de 
TeraWatts ainda poderão ser ainda 
exploradas por outros grupos e ou- 
tros centros de investigação portu- 
gueses da área da Quimica e da 
Biologia. Uma parte relevante do 
tempo de utilização do sistema será 
consagrada a estudos de IDE na 
area da utilização de lasers pulsa- 
dos em projectos industriais. Pre- 
tende-se desta forma, complemen- 
tar a actividade de investigação 
associada a este projecto com acti- 
vidades de Desenvolvimento Expe- 
rimental que resultem naturalmente 
do esforço realizado em torno dos 
projectos de indole fundamental. 


6. CONCLUSÃO 


Neste artigo apresentaram-se os 
conceitos que estão subjacentes 
aos novos mecanismos de acelera- 
ção de partículas carregadas em 
plasmas. As vantagens que estas 
técnicas apresentam foram compa- 
radas com as limitações técnicas e 
financeiras que as actuais experiên- 
cias de altas energias comportam. 
A aceleração de partículas em plas- 
mas supõe a geração de ondas 
electrostáticas de grande amplitude 
e velocidade de fase próximas de c; 
concebem-se actualmente dois me- 
canismos para a criação dessas es- 
truturas em plasmas — beat wave 


e wake field. Falou-se também do 
surfatrão, o esquema proposto para 
evitar o progressivo desfasamento 
entre a partícula carregada e a 
onda electrostática. Finalmente, 
descreveu-se a tecnologia actual de 
lasers de grande potência que per- 
mite a realização destas experiên- 
cias e apresentou-se sumariamente 
o sistema laser a instalar proxima- 
mente no IST, assim como as ca- 
racterísticas genéricas da activida- 
des a desenvolver no âmbito deste 
projecto. 
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Diagrama do Sistema CPA - Tera Watt de Lisboa 
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Impulso Laser 
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Ndvidro 


Figura 2 - Diagrama do sistema CPA-TeraWatt de Lisboa 


RESUMO 


As disciplinas Termodinâmica | 
e Il, leccionadas respectiva- 
mente no 1.º e 2.º semestres 
do 3.º ano da licenciatura em 
Engenharia Mecânica, resulta- 
ram da disciplina Termodina- 
mica, anual, leccionada no 
mesmo 3.º ano até à reforma 
curricular que teve lugar no ini- 
cio da década de 70. Elas 
constituem, no seu conjunto, 
um espaço de aprendizagem 
onde os alunos tomam con- 
tacto com a Termodinâmica 
Macroscopica e as suas princi- 
pais aplicações no domínio da 
Engenharia Mecânica. Note-se 
que os mesmos alunos já 
aprenderam Termodinâmica 
como uma das disciplinas lec- 
cionadas pelo Departamento 
de Fisica, que no entanto é 
tradicionalmente ensinada de 
um ponto de vista microscó- 
pico. 

O presente artigo descreve 
o conteúdo programático da- 
quelas disciplinas, e reflecte 
a experiência do autor como 
responsável e docente das 
mesmas durante vários anos no 
Departamento de Engenharia 
Mecânica do Instituto Superior 
Técnico. Além do programa, 
que se justifica resumidamente, 
apresenta-se uma lista de finali- 
dades e objectivos das discipli- 
nas, e saliegntam-se algumas 
vantagens de um ensino orien- 
tado por objectivos. 

Como os programas das duas 
disciplinas são fortemente in- 
terdependentes, quando não 
der lugar a ambiguidade, ou 
quando forem tratados pontos 
genéricos, referir-se-ã apenas 
Termodinâmica. 


ABSTRACT 
Thermodynamies | and Il are 


one-semester courses, follo- 
wed by all Mechanical Engine- 
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1. INTRODUÇÃO 


Os livros de texto de Termodinã- 
mica para cursos de engenharia, de 
que as referências [1] a [7] consti- 
tuem exemplos típicos, começam 
geralmente, após a exposição de 
alguns conceitos básicos, por enun- 
ciar a 1.º lei da Termodinâmica apli- 
cada a sistemas fechados, pas- 
sando depois à apresentação da 
mesma lei para sistemas abertos. 

A 2.º lei surge posteriormente, e 
com frequência formulada em ter- 
mos do funcionamento das máqui- 
nas térmicas (enunciados de Clau- 
sius e de Kelvin-Planck). 


Verifica-se assim que os conceitos 
de sistema fechado e sistema 
aberto surgem separadamente, e 
de tal modo que problemas envol- 
vendo um e outro tipo de sistemas 
tendem frequentemente a ser consi- 
derados pelos alunos como comple- 
tamente distintos. 


A referência [8] apresenta uma for- 
mulação que difere consideravel- 
mente das anteriores, evitando os 
inconvenientes mencionados, e foi 
assim usada como modelo de parte 
da estrutura do programa que se 
expõe neste artigo. 


Como o tema base da Termodiná- 
mica é a energia e as trocas de 
energia, escolheu-se a noção de 
balanço como a noção central. 


introduzida em relação a uma variá- 
vel extensiva genérica, apresentam- 
se em seguida alguns exemplos 
para mostrar que existem proprie- 
dades conservadas e não conser- 
vadas. Concretizando em relação à 


massa e à energia, são obtidas as 
respectivas equações de balanço 
aplicadas ao caso mais geral (sis- 
tema aberto em regime variável). 
Fazendo as simplificações apropria- 
das nesta equação geral de ba- 
lanço de energia, são deduzidas as 
equações da primeira lei para um 
sistema aberto em regime estacio- 
nário e para um sistema fechado. 
O sistema fechado surge assim 
como um caso particular do sistema 
aberto e não como algo totalmente 
distinto. 


Outro aspecto que se salienta neste 
ponto é o facto de as equações de 
balanço de massa e energia (bem 
como a de entropia, introduzida 
mais tarde) serem apresentadas 
|) referidas a um intervalo de tempo 
e ii) numa base instantânea (em 
termos de taxas de variação, cau- 
dais e potências). 

Esta segunda representação afasta- 
-Se um pouco da Termodinâmica 
Clássica, que apenas trata de esta- 
dos de equilíbrio, onde o tempo não 
intervém como variável; no entanto, 
a sua utilidade nas aplicações justi- 
fica a sua inclusão, e o princípio do 
estado local serve de fundamento a 
essa inclusão. 


A localização desta disciplina no 
curso justifica algumas considera- 
ções adicionais. O 3.º ano é aquele 
em que se pode dizer que os alu- 
nos têm o primeiro contacto com o 
Departamento de Engenharia Me- 
cânica, visto a maior parte das dis- 
ciplinas dos dois primeiros anos 
serem leccionadas pelos departa- 
mentos de Matemática e Física. 
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É assim uma disciplina que, além 
das suas finalidades específicas, 
deve realizar (embora não exclu- 
sivamente) as seguintes tarefas adi- 
cionais: 


e ser um factor de motivação em 
relação ao mundo da engenharia. 
Com este fim, nalgumas aulas 
teóricas são projectados slides ou 
apresentadas transparências com 
instalações e equipamento di- 
verso; 

e ajudar a familiarizar o aluno com 
unidades, ordens de grandeza, 
valores típicos, através da resolu- 
ção de problemas; 

e contribuir para a introdução à 
prática do trabalho experimental, 
mediante a realização de ensaios 
laboratoriais; 


e melhorar a sua capacidade de 
comunicação escrita, através da 
realização de relatórios dos en- 
saios laboratoriais e outros traba- 
lhos, devolvidos após correcção. 


2. PROGRAMA 


Esta secção descreve de forma de- 
talhada o conteúdo das duas disci- 
plinas. Pode ver-se que no primeiro 
semestre se introduz a parte mais 
fundamental — capitulos 1 a 5 — não 
deixando de estabelecer a ponte 
para as aplicações — capítulos 6 a 
13. No segundo semestre, o capi- 
tulo inicial é uma extensão de ma- 
téria já dada no primeiro, enquanto 
nos três restantes capítulos se tra- 
tam aplicações importantes para a 
engenharia mecânica, 


TERMODINÂMICA | 


Capítulo 1 Introdução 


Conceitos básicos. Sistema. Vizi- 
nhança. Propriedades. Hipótese do 
contínuo. Estado. Fase. Equilíbrio. 
Energia. Trabalho. Calor. Estados 
de equilibrio. Sua caracterização. 
Variáveis intensivas e extensivas. 
Equações de estado. 

Estados de não equilibrio. Princípio 
do estado local. 


Capítulo 2 Conservação de Massa e 
Energia 


Equações de balanço. Propriedades 
conservadas e não conservadas. 
Exemplos do quotidiano. Formula- 
ção em termos de um intervalo de 
tempo ou numa base instantânea. 
Equação de balanço para uma pro- 
priedade extensiva genérica. Meto- 
dologia e convenções. 


Equação de balanço de massa. 
Conservação da massa. Casos par- 
ticulares. 

Equação de balanço de energia. 
Conservação de energia. Mecanis- 
mos de transferência de energia. 
Energia associada à massa entrada 
ou saída. Calor. Trabalho. Trabalho 
ao veio. Trabalho de expansão. 
Trabalho associado à transferência 
de massa. Entalpia. Trabalho eléc- 
trico. Casos particulares. 


Capítulo 3 Propriedades das Substân- 
cias 


Algumas propriedades adicionais. 
Energia de Helmholtz, energia de 
Gibbs, entropia. Variáveis naturais. 
Energia interna e propriedades deri- 
vadas em função das suas variaveis 
naturais. Mnemônica do quadrado. 
Relações de Maxwell. Equações de 
estado. Sua determinação. Coefi- 
cientes termodinâmicos: coefi- 
cientes de expansão e de compres- 
sibilidade e calores específicos. 
Relação entre os coeficientes. Jaco- 
bianos e derivadas parciais: método 
das áreas. Propriedades e exem- 
plos de aplicação. 

O modelo de gás perfeito. Equação 
dos gases perfeitos. Constante dos 
gases perfeitos. Algumas implica- 
ções deste modelo. Processos em 
gases perfeitos. Sua representação 
em diversos diagramas termodinã- 
micos. 

Outras equações de estado. Equa- 
ção de Van-der-Waals. Breve refe- 
rência às equações de Bertholet, 
Beattie-Bridgeman, Dieterici, Bene- 
dict-Webb-Rubin, Redlich e Kwong. 
Equação de virial. Factor de com- 
pressibilidade. Lei dos estados 
correspondentes. Diagrama de com- 
pressibilidade generalizado. O mo- 


ering 3rd year students. Consi- 
dered together, they form a 
learning environment where 
students are exposed to Ma- 
croscopic Thermodynamics 
and its main applications in the 
field of Mechanical Enginee- 
ring. The same students have 
already followed a Thermody- 
namics course given by the 
Physics Department, but this is 
usually taught from a micros- 
copic point of view. 

The present paper describes 
the course contents, and re- 
flects the author's experience 
coordinating and teaching 
these courses in the Mechani- 
cal Engineering Department at 
Instituto Superior Técnico for 
several vears. Besides the 
contents, which are briefly jus- 
tified, a list of both courses 
aims and objectives is presen- 
ted, and some advantages of 
an objective oriented approach 
are emphasized. 


delo de líquido incompressível. Al- 
gumas implicações. 

Mudança de fase. Gases e vapo- 
res. Terminologia. Vapor de água. 
Tabelas. Determinação das proprie- 
dades na zona a duas fases. Dia- 
grama temperatura-entropia. Dia- 
grama de Mollier. 


Capítulo 4 Equação de Balanço de En- 
tropia 


Introdução. Processos naturais. En- 
tropia. Equação de balanço de en- 
tropia. Geração de entropia. Casos 
particulares. Balanço de entropia e 
reversibilidade. 

Reversibilidade e irreversibilidade 
em diversos processos. Trabalho 
fornecido e recebido. Rendimentos 
isentrópicos: turbinas, compresso- 
res, tubeiras e difusores. Trabalho 
ao veio reversível. Laminagem. Co- 
eficiente de Joule-Thomson. 


Capítulo 5 Estabilidade e Equilibrio 


Princípio do máximo de entropia e 
do mínimo de energia. Critérios de 
equilíbrio baseados em diferentes 
constrangimentos. Estabilidade in- 
trínseca de sistemas de um só 
componente. Principio de Le Chate- 
lier. Energia de Gibbs. Superfície 
g,P,T, superfície P,v,T e diagramas 
P,T. À água como substância anó- 
mala: algumas consequências. 
Estados metaestáveis. Analogias. 
Equação de Clapeyron. Caso parti- 
cular da transição líquido-vapor. 
Mudança de fase em sistemas com 
mais de um componente. Regra 
das fases de Gibbs. 


Capítulo 6 Ciclos Termodinâmicos 


Motores térmicos. Calor e trabalho 
trocados. Ciclo de Carnot. Rendi- 
mento. Enunciado de Kelvin-Planck 
da 2.º lei. 

Máquinas frigoríficas. Frigoríficos e 
bombas de calor. Enunciado de 
Clausius da 2.º lei. 


Máquinas de combustão interna. 
Aplicação da análise feita para as 
máquinas térmicas às máquinas de 
combustão interna. 


Capítulo 7 Ciclo de Turbina de Gás 


Introdução. Instalações em circuito 
fechado e circuito aberto: vanta- 
gens e desvantagens. 


Análise do ciclo simples. Rendi- 
mento e trabalho útil do ciclo ideal. 
Sua variação com a razão de pres- 
sões. Efeito das irreversibilidades. 
Turbina de gás com regenerador. 
Regenerador ideal e não ideal: efi- 
ciência do regenerador. Estudos 
paramétricos: exemplos. 

Principais caracteristicas das cen- 
trais de turbina de gás existentes 
em Portugal. 


Capítulo 8 Ciclos de Vapor de Produ- 
ção de Potência Mecânica 


Introdução. Distinção entre o ciclo 
de vapor e o ciclo de gás. Esquema 
básico de uma central térmica. Li- 
mitações ao valor do rendimento, 
Parâmetros de caracterização dos 
ciclos de vapor. Rendimentos do ci- 
clo ideal e do ciclo real. Razão de 
rendimentos. Razão de trabalhos. 
Consumo específico de vapor e 
consumo específico de combustível. 
Ciclo de Carnot. Rendimento. Ciclo 
de Rankine. Ciclo de Rankine com 
sobreaquecimento e reaqueci- 
mento. Características do fluido 
ideal para ciclos de vapor. 


Centrais a vapor existentes em Por- 
tugal. Características principais. 


Capitulo 9 Ciclos Frigoríficos 


Introdução. Ciclo de Carnot inverso: 
frigorífico e bomba de calor. Eficiên- 
cias. Efeito refrigerante. Ciclo 
frigorífico de Carnot para um fluido 
condensável. Ciclo de compressão 
de vapor. Diagrama P,h. Conse- 
quências do funcionamento irrever- 
sivel do compressor. Tabelas de 
fluidos frigorigéneos. Propriedades 
desejáveis para estes fluidos. 


Capítulo 10 Outros Ciclos 


Ciclos combinados: razões para a 
sua utilização. Ciclo combinado 
mercúrio-vapor de água. Ciclo com- 
binado gás-vapor. Motores de com- 
bustão interna. Ciclos de ar padrão. 
Ciclos de Otto, Diesel e misto. Defi- 
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nições e rendimentos. Diagramas 
Pv para motores reais. Ciclo de Stir- 
ling: ciclo ideal e realização prática. 


Capítulo 11 Cogeração 


Introdução. Definição. Vantagens 
da cogeração, Turbina de contra- 
pressão. Turbina de extracção-con- 
densação. Unidades de pequena 
potência: o TOTEM — característi- 
cas principais. 


Capitulo 12 Recuperação de Energia 


Análise do uso da energia. Pontos 
principais a considerar. Energia per- 
dida. Origens possíveis dessas per- 
das. Análise custo-benefício. Perdas 
de energia. Qualidade e quantidade. 
Sistemas de recuperação de calor. 
Métodos de recuperação. 


Capitulo 13 Auditorias de Energia 


Esquema geral da auditoria. Regis- 
tos de consumo. Manutenção. Pes- 
soal. Investimentos. Diferentes for- 
mas de apresentar a análise 
energética. Histogramas. Diagramas 
circulares. Diagrama de Sankey. 
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Capítulo 1 Ciclos Avançados 


Ciclo de turbina de gás. Reaqueci- 
mento. Arrefecimento intermédio. 
Influência no trabalho útil. Sua utili- 
zação em conjunto com regenera- 
ção. Optimização do rendimento. 
Ciclos de vapor. Ciclos com rege- 
neração. Subtiragens de vapor. 
Pré-aquecedores abertos e fecha- 
dos. Economizador e pré-aquece- 
dor de ar. Ciclos de vapor em cen- 
trais nucleares. 


Ciclos frigoríficos. Compressão em 
dois andares; tanque de flash. Ci- 
clos de compressão de vapor em 
cascata. Aplicações. Liquefacção 
de gases: processo de Linde. Ciclo 
frigorífico de gás. refrigerador a 
agua. Frigorífico de absorção. 

Breve referência ao funcionamento 
de instalações a carga parcial. 
Curvas de funcionamento de alguns 


componentes. Consequências do 
funcionamento da instalação em 
condições diferentes da carga no- 
minal. Exemplos. 


Capítulo 2 Misturas Gasosas Ideais. 
Psicrometria 


Misturas e misturas ideais. Compo- 
sição de uma mistura. Fracções 
mássica e molar. Leis de Dalton e 
Amagat-Leduc. Misturas de gases 
perfeitos. O ar húmido. Humidade 
específica, humidade relativa e 
ponto de orvalho. Entalpia do ar hú- 
mido. Temperatura húmida e seca. 
Diagrama psicrométrico. 

Processos com ar húmido. Aqueci- 
mento, arrefecimento, humidifica- 
ção, desumidificação, mistura de 
caudais. Balanços mássicos e ener- 
gético. Equipamentos. Factor de 
contacto. Condicionamento de ar. 
Conforto térmico. Torres de arrefe- 
cimento. 


Capítulo 3 Introdução à Combustão 


Importância da combustão. Com- 
bustíveis: sua caracterização. Bre- 
ves noções sobre a tecnologia de 
queima. 

Reacções quimicas. Estequiome- 
tria. Equações químicas. Análise 
dos produtos da combustão. 


Análise energética da combustão. 
Combustão a volume constante e 
a pressão constante. Entalpia de 
combustão e de formação. Poder 
calorífico. Temperatura de com- 
bustão adiabática: seu cálculo. 


Equilíbrio. Constante de equilíbrio. 
Princípio de Le Chatelier. Dissocia- 
ção. Composição de equilíbrio de 
uma mistura com dissociação. Tem- 
peratura de combustão adiabática 
quando ocorre dissociação. 


Capítulo 4 Exergia 


Introdução. Conceito de exergia. 
Ambiente, equilíbrio, estado morto. 
Exergia associada à transferência 
de energia sob a forma de trabalho, 
de calor e a um escoamento em re- 
gime estacionário. Componentes da 
exergia: cinética, potencial, física e 
quimica. 


Análise exergética de sistemas fe- 
chados. Exergia estática. Equação 
de balanço de exergia. Casos 
particulares. Aplicação da equação 
de balanço de exergia a sistemas 
abertos em regime estacionário. 
Taxa de irreversibilidade. Rendi- 
mento exergético. Representação 
gráfica do balanço de exergia. 


Exergia química. Substâncias de re- 
ferência. Exergia química de uma 
mistura. Análise exergética de 
processos simples. Expansão, lami- 
nagem, compressão, combustão. 
Formas de diminuir a irreversibili- 
dade. Exergia química de combusti- 
veis industriais. 


3. FINALIDADES E OBJECTIVOS 


O ensino orientado por objectivos é 
tratado de forma exaustiva em 
Bloom [9]. A ideia central do mé- 
todo é que qualquer programa de 
instrução tem como finalidade aju- 
dar os estudantes a mudar o seu 
comportamento; após a instrução, O 
estudante deve ser capaz de fazer 
alguma coisa que antes não era ca- 
paz. 


Uma vez que os objectivos sejam 
definidos claramente, podem tornar- 
-se modelos ou planos que ajudam 
a formar e a guiar os processos de 
instrução e avaliação. 

A descrição dos objectivos deve ser 
feita em termos operacionais, e as 
palavras usadas não devem ser do 
tipo conhecer / compreender / pos- 
suir espírito crítico / apreciar... que 
se caracterizam por uma certa 
imprecisão e são aptidões que não 
é possível medir, mas sim do tipo 
declarar / dizer pelas próprias pala- 
vras / calcular / definir... 

A presente secção expõe as finali- 
dades e objectivos das disciplinas. 
As finalidades indicam as principais 
linhas de força da disciplina, de 
uma forma necessariamente con- 
densada. Os objectivos enunciam 
capacidades que se pretende que o 
aluno possua após o estudo de 
cada ponto da matéria; são assim 
formulados com o máximo por- 
menor e sempre em termos de um 
comportamento que o aluno deve 
mostrar. 


A lista apresentada pode ter duas 
utilizações: 


e pode servir ao aluno para uma 
auto-avaliação antes de compare- 
cer a um exame; 

* pode ser um auxiliar na normali- 
zação das classificações: a partir 
dela é possível seleccionar um 
subconjunto de objectivos que 
constitua o mínimo exigivel para 
aprovação na disciplina. 


Se cada disciplina leccionada no 
IST fizesse este trabalho, e garan- 
tisse que cada aluno passado tinha 
satisfeito os objectivos mínimos, 
tornando ao mesmo tempo con- 
hecidos esses objectivos mínimos 
as disciplinas a juzante, o problema 
das precedências seria desdramati- 
zado e poderia ser discutido numa 
base racional e objectiva. 


Finalidades 


Familiarizar o aluno com os concei- 
tos de energia e de entropia e as 
formas como estas grandezas po- 
dem ser transferidas e armazena- 
das; 


Torná-lo capaz de aplicar com faci- 
lidade o formalismo termodinâmico, 
bem como de utilizar a informação 
sobre propriedades termodinâmicas 
(tabelas ou diagramas). 


Formular e aplicar balanços de 
massa, de energia e de entropia. 
Apresentar exemplos diversos de 
sistemas termodinâmicos, forne- 
cendo uma metodologia geral de 
análise que permita ao aluno o es- 
tudo de qualquer processo ou sis- 
tema termodinâmico. 

Introduzir o conceito de exergia, for- 
mular a equação de balanço de 
exergia e aplicá-la na análise de 
sistemas energéticos. 


Objectivos 


Por uma questão de sistematiza- 
ção, os objectivos que se seguem 
foram distribuidos pelos dois se- 
mestres lectivos e, dentro de cada 
um, pelos capítulos em que está di- 
vidida a matéria. Assim, depois de 
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estudar cada capítulo, o aluno deve 
ser capaz de: 


TERMODINÂMICA | 


Capítulo 1 Introdução 


Definir sistema, vizinhança, proprie- 
dade, estado e fase. 

Caracterizar o estado de equilíbrio 
de um sistema termodinâmico. 
Definir variável intensiva e exten- 
siva e apresentar exemplos. 

Definir variável de estado e equa- 
ção de estado. Exemplificar. 


Enunciar o princípio do estado local 
e justificar a sua importância. 


Capítulo 2 Conservação de Massa e 
Energia 


Formular uma equação de balanço 
de uma propriedade extensiva gené- 
rica, quer em termos da variação 
instantânea quer da variação num 
intervalo de tempo finito, e caracteri- 
zar os diversos termos da equação. 
Distinguir entre propriedade conser- 
vada e propriedade não conser- 
vada. Dar exemplos. 


Reconhecer a diferença entre a 
equação de balanço de uma e outra. 


Escrever a equação de balanço de 
massa para um sistema aberto em 
regime não estacionário. 
Simplificar a equação para tratar o 
caso de regime estacionário e sis- 
tema fechado. 


Reconhecer as componentes da 
energia total de um sistema termo- 
dinâmico (interna, cinética e poten- 
cial). 

Descrever os diversos mecanismos 
de transferência de energia entre 
um sistema e a sua vizinhança: 
fluxo de energia associado a um 
fluxo de massa, calor, trabalho 
(trabalho ao veio, trabalho de ex- 
pansão, trabalho associado à trans- 
ferência de massa, trabalho eléc- 
trico). 

Escrever a equação de balanço de 
energia para o caso geral de um 
sistema aberto em regime não es- 
tacionário e identificar todos os ter- 
mos. 


Simplificar a equação para tratar o 
caso de sistema adiabático, sistema 
aberto em regime estacionário e 
sistema fechado. 


Capítulo 3 Propriedades das Substân- 
cias 


Caracterizar uma variável natural; 
exprimir as diferenciais da energia 
interna U, entalpia H, energia de 
Helmholtz F e energia de Gibbs G 
em função das diferenciais das 
suas variáveis naturais. 


Obter as relações de Maxwell a par- 
tir das expressões do ponto anterior. 


Utilizar a mnemónica do quadrado 
para obter as expressões dos dois 
pontos anteriores. 

Definir os coeficientes de expansão 
isobárica e de compressibilidade 
isotérmica, e os calores específicos 
a pressão e a volume constante; 
dada a equação que relaciona es- 
tes quatro coeficientes, verificar que 
se trata de uma identidade. 


Interpretar o jacobiano de uma 
transformação entre dois espaços 
bi-dimensionais como um quociente 
de áreas elementares; utilizar o mé- 
todo das áreas para exprimir uma 
derivada parcial em que interve- 
nham quaisquer três propriedades 
em função dos coeficientes ter- 
modinâmicos. 

Exprimir as variações de energia 
interna e entalpia dU e dH em fun- 
ção das variações de quaisquer 
duas das três variáveis P, Ve T. 
Escrever a equação dos gases per- 
feitos (em termos molares ou más- 
sicos e com o volume total ou vo- 
lume específico) e utilizá-la. 

Obter a constante dos gases perfei- 
tos R para qualquer gás a partir da 
massa molar e da constante univer- 
sal RO e conhecer o valor desta. 
Definir processos isócoro, isobárico, 
adiabático, isotérmico e representá- 
los nos diagramas P,T; Pv e Ty. 
Mostrar que a energia interna e a 
entalpia de um gás perfeito são 
apenas função da temperatura. 
Determinar a variação da entropia 
de um gás perfeito em função das 
variações de quaisquer duas das 
variáveis Pv e T. 


Escrever as equações de uma adi- 
abática e de uma politrópica e co- 
nhecer as condições de aplicação. 
Resolver problemas envolvendo 
processos simples com gases per- 
feitos. 

Escrever a equação de Van der 
Waals, identificar o que a distingue 
da equação dos gases perfeitos e 
inferir a forma das isotérmicas num 
diagrama Pv. 


Definir factor de compressibilidade 
e enunciar a lei dos estados corres- 
pondentes. 

Caracterizar o modelo de liquido in- 
compressivel. Obter expressões 
para a variação de energia interna 
e entalpia de um líquido incompres- 
sivel. 

Caracterizar uma mudança de fase 
(sólido-líquido, liquido-vapor, sólido- 
-vapor) e definir gás e vapor. 

Estar familiarizado com a terminolo- 
gia da mudança de fase. 


Definir calor latente, pressão de 
saturação, temperatura de satura- 
ção e grau de sobreaquecimento. 


Usar com facilidade as tabelas de 
vapor de água e de fluidos frigori- 
géêneos. 

Definir título do vapor e usá-lo para 
calcular propriedades na zona a 
duas fases. 


Esboçar nos diagramas T,s e h,s as 
isolinhas mais importantes para um 
fluido condensável. 


Capítulo 4 Equação de Balanço de En- 
tropia 


Definir entropia e escrever a sua 
equação de balanço, em termos 
instantâneos e para um intervalo de 
tempo finito, e identificar todos os 
termos. 


Definir processo reversível e identi- 
ficar as condições em que um pro- 
cesso pode ser considerado como 
tal. 


Mostrar que o trabalho produzido 
por um sistema é máximo para um 
processo reversível. 

Mostrar que o trabalho recebido por 
um sistema é minimo para um pro- 
cesso reversível. 


Definir rendimento isentrópico de 
uma turbina e de um compressor 
(ou bomba); mostrar o significado 
deste parâmetro para aferir o com- 
portamento de um sistema real, 


Caracterizar uma tubeira e um difu- 
sor e explicar a sua utilidade. Apli- 
car-lhes as equações de balanço 
mássico e energético; definir os 
respectivos rendimentos isen- 
trópicos e mostrar o seu significado. 


Definir laminagem de um fluido e 
aplicar a uma válvula de laminagem 
a equação de balanço de energia. 
Esboçar as isentálpicas e a curva 
de inversão num diagrama P,T; de- 
finir o coeficiente de Joule-Thomson 
e mostrar a sua utilidade. 


Capítulo 5 Estabilidade e Equilíbrio 


Enunciar o princípio do máximo de 
entropia e o princípio do mínimo de 
energia. 

Deduzir a condição de equilíbrio 
para um sistema fechado com Te V 
constantes e com T e P constantes. 
Mostrar que o critério de estabi- 
lidade para um sistema de um só 
componente é expresso por 

| OP) 

usa E ) <() 

Enunciar o princípio de Le Chatelier 
e dar exemplos da sua aplicação. 
Esboçar as linhas de coexistência 
de fases no diagrama P,T; distinguir 
o caso da água do da generalidade 
das substâncias. Enunciar conse- 


quências do comportamento anó- 
malo da água. 


Definir Estado Metaestável e Exem- 
plificar Situações em Que Tais Es- 
tados Podem Ocorrer. 

Deduzir a equação de Clapeyron; 


aplicá-la à transição líquido-vapor 


para obter uma expressão aproxi- 
mada da pressão do vapor em fun- 
ção da temperatura. 

Enunciar a regra das fases de 


Gibbs e exemplificar a sua aplica- 
ção. 


Capitulo 6 Ciclos Termodinâmicos 


Definir ciclo. Calcular o trabalho e o 
calor postos em jogo num ciclo e 
mostrar a sua correspondência com 
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as áreas limitadas pelos ciclos nos 
diagramas Pve Ts. 

Definir motor térmico; definir fonte 
quente e fonte fria e dar exemplos 
de concretização destes dois con- 
ceitos. 


Definir ciclo de Carnot, explicar a 
sua utilidade, representá-lo nos dia- 
gramas P,v, T;se Tv. 


Definir rendimento e exprimi-lo em 
termos das temperaturas da fonte 
quente e fonte fria, para um motor 
reversível e tendo em conta as irre- 
versibilidades. 


Conhecer o enunciado de Kelvin- 
Planck da 2.º lei da Termodinâmica 
e demonstrá-lo a partir da equação 
de balanço de entropia. 

Definir máquina frigorifica; distinguir 
entre frigorífico e bomba de calor; 
definir eficiência em ambos os casos. 


Conhecer o enunciado de Clausius 
da 2.º lei e demonstrá-lo a partir da 
equação de balanço de entropia. 
Provar a equivalência dos enuncia- 
dos de Clausius e de Kelvin-Planck. 
Listar as hipóteses simplificativas 
que permitem analisar máquinas de 
combustão interna como motores 
termicos. 


Capítulo 7 Ciclo de Turbina de Gás 


Desenhar esquematicamente uma 
instalação de turbina de gás tanto 
em circuito aberto como em circuito 
fechado, identificando os compo- 
nentes principais. Reconhecer van- 
tagens e desvantagens de um e ou- 
tro tipo de instalação. 

Representar o ciclo nos diagramas 
n,s; [;se Pu. 


Obter a expressão do rendimento 
do ciclo ideal em função da razão 
de pressões. 

Obter uma expressão para o traba- 
lho útil do ciclo e identificar o valor 
da razão de pressões para o qual 
o trabalho util é máximo. 

Indicar as principais consequências 
do funcionamento irreversível de 
um compressor ou turbina. 


Obter expressões para o rendi- 
mento e para o trabalho útil no 
caso do ciclo irreversível, identificar 
as diferenças em relação ao ciclo 
reversivel. 


Explicar a função do regenerador 
numa instalação de turbina de gás 
e conhecer as condições de utiliza- 
ção. Definir eficiência do regenera- 
dor e mostrar a influência do seu 
valor no rendimento do ciclo com 
regeneração. 

Resolver problemas numéricos 
usando a teoria dos pontos ante- 
rRores. 


Capítulo 8 Ciclos de Vapor de Produ- 
ção de Potência Mecânica 


Distinguir entre ciclo de gás e ciclo 
de vapor. 

Definir os seguintes parâmetros uti- 
lizados na caracterização dos ciclos 
de vapor: rendimento do ciclo ideal, 
rendimento do ciclo real, razão de 
rendimentos, razao de trabalhos, 
consumo especifico de vapor e 
consumo específico (de combusti- 
vel) da central. 


Caracterizar o ciclo de Carnot fun- 
cionando com um fluido condensá- 
vel e explicar porque não é possível 
de realizar na prática. 

Caracterizar o ciclo de Rankine e 
mostrar as suas vantagens em rela- 
ção ao ciclo de Carnot funcionando 
entre as mesmas temperaturas. 
Definir sobreaquecimento e reaque- 
cimento e indicar as respectivas 
vantagens. 


Resolver problemas numéricos utili- 
zando as tabelas de vapor de água. 


Capítulo 9 Ciclos Frigoríficos 


Descrever o funcionamento de um 
frigorífico e de uma bomba de calor 
baseados no ciclo frigorífico de Car- 
not; definir e obter expressões para 
as respectivas eficiências. Definir 
efeito frigorífico. 


Descrever um ciclo frigorífico de 
Carnot para um fluido condensável 
e representá-lo no diagrama Ts. 
Enunciar as principais razões pelas 
quais este ciclo não é usado. 


Descrever o ciclo de compressão 
de vapor, representá-lo no dia- 
grama T,s e saber as diferenças 
entre este ciclo e o ciclo frigorífico 
de Carnot. 


Esboçar um diagrama P,h com a lo- 
calização das principais isolinhas e 


representar neste diagrama um ci- 
clo de compressão de vapor. 


Resolver problemas de ciclos frigo- 
ríficos usando as tabelas ou diagra- 
mas. 


Capítulo 10 Outros Ciclos 


Indicar a principal vantagem da utili- 
zação de ciclos combinados. Des- 
crever as principais vantagens e 
inconvenientes do ciclo combinado 
mercúrio/vapor de agua. Descrever 
as características principais de uma 
instalação combinada turbina de 
gás/turbina de vapor. 

Descrever os ciclos de ar padrão 
Otto, Diesel e misto e obter expres- 
sões para os respectivos rendi- 
mentos. 


Esboçar o diagrama P,v de um mo- 
tor a gasolina e de um motor Die- 
sel; indicar as diferenças em rela- 
ção aos ciclos de ar padrão 
correspondentes. 

Descrever o ciclo de Stirling; repre- 
sentá-lo no diagrama P,v; referir os 
principais problemas que se colo- 
cam à sua realização prática. 


Capítulo 11 Cogeração 


Definir cogeração, explicar o seu in- 
teresse e apresentar casos em que 
se justifique O seu uso. 


Descrever algumas instalações de 
cogeração, tanto de média/alta po- 
tência como de pequena potência. 


Capítulo 12 Recuperação de Energia 


Identificar diversas formas como, 
num processo ou instalação indus- 
trial, a energia pode ser perdida. 


Descrever os métodos de recupe- 
ração de calor mais utilizados. 


Capitulo 13 Auditorias de Energia 


Saber em que consiste uma audito- 
ria de energia, qual a sua finali- 
dade, quais os principais passos 
que a constituem. 

Conhecer as diversas formas de 
apresentar os resultados de uma 
análise energética. 


TERMODINÂMICA II 


Capítulo 1 Ciclos Avançados 


Ciclo de turbina de gás 


Caracterizar o reaquecimento e 
mostrar o seu interesse; explicar a 
sua exequibilidade e os seus limites 
em instalações de circuito aberto. 
Caracterizar a compressão com ar- 
refecimento intermédio. Calcular a 
pressão intermédia que maximiza o 
trabalho útil. Mostrar o interesse de 
usar estas duas modificações em 
conjunto com a regeneração. 


Resolver problemas numéricos de 
aplicação da teoria do ponto ante- 
ror. 


Ciclos de vapor 


Descrever a regeneração e explicar 
as suas vantagens. Discutir as con- 
sequências das subtiragens de va- 
por. Explicar como se faz na prática 
a regeneração, distinguindo entre 
pré-aquecedores abertos e fecha- 
dos e as implicações da utilização 
de uns ou outros. Resolver proble- 
mas de ciclos de vapor, Incluindo 
sobreaquecimento, reaquecimento 
e regeneração. 

Explicar a função do economizador 
e do pré-aquecedor de ar numa 
central térmica. 


Ciclos frigoríficos 


Descrever e representar esquemati- 
camente um ciclo frigorífico com 
compressão em dois andares e uti- 
lização de um tanque de flash; 
idem para uma instalação frigorífica 
com ciclos de compressão de vapor 
em cascata. Resolver problemas 
numéricos em instalações deste 
tipo utilizando as tabelas de fluidos 
frigorigéneos. 


Descrever o processo de Linde de 
liquefacção de gases. 

Descrever o ciclo frigorífico de gás 
e listar as suas vantagens e incon- 
venientes quando comparado com 
o ciclo de compressão de vapor. 
Descrever o ciclo frigorífico a água 
(utilizando uma bomba de ejector) 
Explicar o princípio de funciona- 
mento e descrever sucintamente 
um frigorífico de absorção. 
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Funcionamento a carga parcial 


Explicar a diferença entre o funcio- 
namento de uma instalação à carga 
nominal e a carga parcial, com re- 
curso a curvas de funcionamento 
dos vários componentes. 


Capítulo 2 Misturas Gasosas Ideais. 
Psicrometria 


Definir mistura e mistura ideal. 
Definir fracção mássica, fracção 
molar e fracção volumétrica. 
Enunciar as leis de Dalton e Ama- 
gat-Leduc e dar o seu significado fi- 
Sico. 


Passar de fracções massicas a 
fracções molares e inversamente. 


Mostrar que o vapor de água con- 
tido no ar pode ser considerado um 
gás perfeito. 

Definir humidade específica, humi- 
dade relativa e ponto de orvalho. 


Derivar a expressão que dá a en- 
talpia do ar húmido em função das 
entalpias específicas do ar seco e 
vapor de água e da humidade es- 
pecífica. 

Definir temperatura húmida e tem- 
peratura seca. 

Descrever um termómetro de bolbo 
húmido, seu princípio de funciona- 
mento e diferença em relação ao 
termómetro de bolbo seco. 


Esboçar um diagrama psicromé- 
trico, incluindo nele isolinhas de hu- 
midade relativa, temperatura seca, 
temperatura húmida e volume espe- 
cifico. 


Em relação aos seguintes proces- 
sos com ar húmido 


a) aquecimento 


b) arrefecimento 

c) humidificação 

d) desumidificação 
e) mistura de caudais 


descrever os processos e repre- 
sentá-los no diagrama psicromé- 
trico; 

efectuar os balanços mássicos e 
energético a um volume de controlo 
onde o processo tem lugar; 


descrever os equipamentos onde 
os processos tem lugar, em parti- 
cular baterias de aquecimento e ar- 


refecimento e lavadores de ar, com 
e sem temperatura controlada. 


Definir factor de contacto e conhe- 
cer valores típicos deste parâmetro. 


Calcular as caracteristicas do cau- 
dal de ar húmido resultante da mis- 
tura de dois caudais dados. 


Descrever em que consiste o condi- 
cionamento de ar e quais as pro- 
priedades do ar que interessa con- 
trolar. 


Justificar a necessidade do condi- 
cionamento de ar em termos do ba- 
lanço energético de um organismo 
vivo; definir e caracterizar conforto 
térmico. 


Descrever todos os componentes 
de um condicionador de ar para po- 
der efectuar todas as operações 
que eventualmente sejam necessá- 
rias (com tempo quente ou frio). 


Descrever uma torre de arrefeci- 
mento, explicar o seu funciona- 
mento e utilidade. 


Capítulo 3 Introdução à Combustão 


Definir combustão , fogo e chama, 
Dar exemplos de sistemas e pro- 
cessos em que a combustão seja 
essencial. 

Caracterizar o ar atmosférico atra- 
vés das proporções dos seus prin- 
cipais constituintes. 

Dar exemplos de combustíveis sóli- 
dos. 


Definir análise elementar e análise 
imediata de um carvão. 


Dar exemplos de combustíveis li- 
quidos. 

Caracterizar os combustiveis líqui- 
dos mais correntes através do hi- 
drocarboneto predominante. 


Dar exemplos de combustiveis ga- 
sosos. 

Caracterizar a diferença entre cha- 
mas de pré-mistura e chamas de 
difusão , incluindo os aspectos liga- 
dos à segurança. 

Indicar os principais problemas as- 
sociados à queima dos combusti- 
veis liquidos, referindo a importân- 
cia da atomização do combustível. 


Descrever as diversas formas de 
queimar combustível sólido, refe- 


rindo nomeadamente a combustão 
em grelha, sob a forma de combus- 
tível pulverizado e em leito fluidi- 
zado. 


Definir estequiometria, reagentes e 
produtos de uma reacção. 


Definir equação quimica. 
Escrever as equações químicas da 
combustão dos combustíveis mais 


correntes, com oxigénio e ar, numa 
base molar e mássica. 


Definir ar estequiométrico ou teó- 
rico, excesso de ar, mistura pobre, 
mistura rica, razão de equivalência, 
riqueza da mistura, combustão 
neutra, oxidante e redutora. 


Indicar a ordem de grandeza do ex- 
cesso de ar para os processos de 
combustão mais correntes. 


Definir análise húmida e seca dos 
produtos da combustão. 

Escrever a equação de balanço 
energético para um sistema fe- 
chado a volume constante, sistema 
aberto a pressão constante e sis- 
tema fechado a pressão constante, 
dentro dos quais tem lugar uma 
reacção de combustão. Particulari- 
zar para os casos extremos em 
que: 1) o processo é adiabático; 
2) a temperatura final dos produtos 
da combustão é igual à tempera- 
tura inicial dos reagentes e repre- 
sentá-los, bem como um caso 
intermédio, num diagrama energia- 
“temperatura. 

Definir temperatura de combustão 
adiabática. 

Definir entalpia (energia interna) de 
combustão e entalpia (energia in- 
terna) de combustão padrão. 

Dada a entalpia (energia interna) de 
combustão padrão, calcular a en- 
talpia (energia interna) de combus- 
tão a qualquer outra temperatura. 
Definir entalpia de formação e en- 
talpia de formação padrão. 

Definir poder calorífico de um com- 
bustivel; distinguir entre poder calo- 
rífico a pressão constante e a vo- 
lume constante e entre poder 
calorífico superior e inferior. 
Calcular a temperatura de combus- 
tão adiabática, usando 1) tabela de 
calores específicos em função da 


temperatura; 2) tabela de entalpias 
em função da temperatura. 


Definir constante de equilibrio; ex- 
primi-la em termos de pressões 
parciais e fracções molares e rela- 
cionar entre si estas duas expres- 
sões. 

Enunciar o princípio de Le Chate- 
lier; utilizá-lo para explicar a influên- 
cia no equilíbrio de 1) variação da 
concentração de um dos constituin- 
tes; 2) variação do volume total (ou 
pressão total) do sistema; 3) varia- 
ção da temperatura do sistema. 
Definir dissociação e descrever qual 
o seu efeito em combustão . 


Calcular a composição de uma mis- 
tura de produtos típicos de uma 
reacção de combustão, dada a tem- 
peratura e uma tabela de constantes 
de equilibrio. 

Calcular a temperatura de combus- 
tão adiabática, considerando a dis- 
sociação dos produtos. 


Capítulo 4 Exergia 


Caracterizar exergia em relação a 
energia. Definir ambiente, equilíbrio 
restrito com o ambiente, estado 
morto. 


Definir a exergia associada ao tra- 
balho, ao calor e a um escoamento. 
Identificar no último caso as parce- 
las cinética, potencial, fisica e qui- 
mica da exergia. 

Mostrar que a exergia física tem 
uma componente térmica e uma 
componente de pressão. 


Escrever a expressão da exergia 
contida num sistema fechado. 


Escrever a equação de balanço de 
exergia aplicada a um sistema 
aberto em regime não estacionário, 
identificar todos os termos e aplicá- 
la a problemas concretos. 

Calcular a taxa de irreversibilidade 
em processos simples. 


Definir rendimento exergético e cal- 
cular o seu valor para uma dada 
instalação. Explicar o que o distin- 
gue do rendimento energético. 


Conhecer e aplicar os diversos mé- 


todos de representação gráfica do 


balanço de exergia. 
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Definir exergia química, identificar 
as substâncias de referência, calcu- 
lar a exergia química de uma mis- 
tura. 

Efectuar balanços de exergia a di- 
versos equipamentos, nomeada- 
mente onde ocorram processos de 
expansão, compressão, laminagem 
ou de combustão. 

Identificar formas de diminuir a irre- 
versibilidade dos processos. 
Calcular a exergia química de com- 
bustíveis industriais. 


4. JUSTIFICAÇÃO DO 
CONTEÚDO 


Na presente secção incluem-se al- 
guns comentários sobre o conteúdo 
do programa. As razões da inclusão 
de um determinado assunto deri- 
vam ou da coerência interna do 
programa, ou de uma necessidade 
(detectada ou antecipada) na prá- 
tica profissional de um engenheiro 
mecânico. 


TERMODINÂMICA | 


O Capítulo 1 destina-se a introdu- 
ção de alguns conceitos básicos e 
de um conjunto de termos que vão 
constituir o vocabulário da disci- 
plina. Esta introdução é relati- 
vamente rápida, uma vez que se 
presume que os conceitos são já 
conhecidos da disciplina Termodi- 
nâmica do 2.º ano (embora se veri- 
fique que, mais vezes do que seria 
de desejar, isso não acontece). Ca- 
racterizam-se os estados de equili- 
brio dos sistemas termodinâmicos, 
apresenta-se o conceito de equa- 
ção de estado e através da introdu- 
ção do princípio do estado local, 
justifica-se a aplicação das relações 
de equilibrio a estados de não 
equilibrio, sendo apontados os limi- 
tes desta aproximação. 

O Capítulo 2 introduz a noção de 
balanço. Apresentadas algumas 
aplicações com exemplos do quoti- 
diano, mostra-se a diferença entre 
propriedades conservadas e não 
conservadas e formula-se a equa- 
ção de balanço para uma proprie- 
dade extensiva genérica. Esta 


equação é em seguida concretizada 
para a massa e energia, definindo 
cuidadosamente os diversos termos 
que entram nestas equações de ba- 
lanço. Estas equações são sempre 
apresentadas na sua forma mais 
geral — aplicadas a um sistema 
aberto em regime não estacionário 
— sendo depois, mediante simplifi- 
cação de certos termos, deduzidos 
diversos casos particulares de 
maior interesse, nomeadamente o 
sistema aberto em regime estacio- 
nário e o sistema fechado. Con- 
forme se mencionou na Introdução, 
as equações de balanço são tam- 
bém apresentadas em duas ver- 
sões: referidas a um intervalo de 
tempo e numa base instantânea. 
Com a finalidade de permitir quanti- 
ficar os problemas a resolver, o Ca- 
pítulo 3 apresenta a metodologia 
que permite relacionar as proprie- 
dades termodinâmicas entre si, in- 
troduz os coeficientes termodinâmi- 
cos e apresenta os casos mais 
típicos de equação de estado: a 
equação dos gases perfeitos (que 
permite resolver problemas em que 
o ar e os chamados gases per- 
manentes são os fluidos de traba- 
lho) e equações de estado expres- 
sas por tabelas ou diagramas 
(discute-se a fenomenologia da mu- 
dança de fase, concretizando para 
os casos do vapor de água e flui- 
dos refrigerantes). Chama-se a 
atenção para as limitações do mo- 
delo de gás perfeito, referindo as 
equações de estado para gases 
não perfeitos, mas dando uma 
atenção especial à equação de Van 
der Waals, tanto pela sua simplici- 
dade como por servir de ponte para 
o fenómeno de mudança de fase. 
Refere-se também neste contexto o 
factor de compressibilidade e a lei 
dos estados correspondentes. 
É ainda apresentado o modelo de 
liquido incompressível. 

O Capítulo 4 define entropia e apre- 
senta a sua equação de balanço; o 
conceito de processo reversível 
surge neste capítulo como o caso 
limite para o qual a geração de en- 
tropia é igual a zero. Associa-se a 
reversibilidade ao máximo (minimo) 
do trabalho fornecido (recebido) por 
um sistema que evolui entre condi- 


ções dadas. Descrevem-se os pro- 
cessos que ocorrem em turbinas, 
compressores, tubeiras, difusores e 
válvulas de laminagem, introdu- 
zindo-se o conceito de rendimento 
isentrópico. Deduz-se também a ex- 
pressão para o trabalho ao veio re- 
versível num sistema aberto em re- 
gime estacionário. 

O Capítulo 5 trata com mais profun- 
didade o conceito de equilíbrio, já 
introduzido no capítulo 1, bem 
como a noção de estabilidade de 
um sistema. Discute-se a esta- 
bilidade intrinseca de sistemas de 
um só componente, chamando-se a 
atenção para o comportamento 
anómalo da água. A mudança de 
fase é apresentada como o exem- 
plo de uma forma que um sistema 
tem de “fugir” à instabilidade. Es- 
tuda-se particularmente a transição 
líquido-vapor, de maior interesse 
em engenharia mecânica. O capí- 
tulo termina com a apresentação da 
regra das fases de Gibbs. 

Uma aplicação muito importante da 
Termodinâmica em engenharia me- 
cânica são as máquinas de funcio- 
namento cíclico. Os capítulos se- 
guintes dedicam-se a este aspecto. 
O conceito de ciclo termodinâmico 
é apresentado no Capítulo 6 como 
O processo que ocorre num sistema 
que passa por uma sucessão de 
estados, regressando ao estado ini- 
cial, enquanto troca calor com duas 
fontes (reservatórios de energia tér- 
mica) a temperaturas diferentes. 
Este nível de generalidade permite 
definir máquina motriz e máquina 
frigorífica. Faz-se a distinção entre 
máquina térmica como definida 
acima e máquina de combustão in- 
terna, e apresentam-se as aproxi- 
mações usuais para analisar estas 
com a metodologia desenvolvida 
para a análise dos ciclos. 

Após esta introdução genérica aos 
ciclos, os capítulos 7 a 10 apresen- 
tam diversos tipos particulares de 
ciclo termodinâmico: o ciclo de tur- 
bina de gás (Capítulo 7), o ciclo de 
vapor de produção de potência me- 
cânica (Capítulo 8) e o ciclo frigori- 
fico (Capítulo 9) são tratados a um 
nível relativamente elementar, e se- 
rão de novo abordados em Termo- 
dinâmica Il. São ainda dados, de 


uma forma mais qualitativa, alguns 
exemplos de ciclos combinados 
(mercúrio/vapor de água e gás/ 
tâágua), os ciclos Otto, Diesel e 
misto e o ciclo de Stirling (Capítulo 
10). Deste último será feita uma de- 
monstração na aula teórica, en- 
quanto os ciclos de turbina de gás, 
vapor e frigorífico serão objecto de 
ensaios laboratoriais. 

O Capítulo 11 é uma última aplica- 
ção de ciclos termodinâmicos: tra- 
tando de cogeração, nele se des- 
crevem unidades de média potência 
(instalações de vapor) e de pe- 
quena potência, baseadas no motor 
de automóvel. Os capítulos 12 e 13 
tratam de aspectos concretos de 
utilização da análise energética. No 
Capítulo 12, são listadas diversas 
origens do desperdício de energia e 
descritos vários metodos de recu- 
peração de calor. Esta problemática 
é retomada de forma mais geral no 
Capítulo 13, onde se descrevem 
com certo pormenor os diversos 
passos que constituem uma audito- 
ria de energia. 


TERMODINÂMICA Il 


O Capítulo 1 é a continuação dos 
capítulos 7, 8 e 9 do semestre an- 
terior. São estudados de forma 
mais detalhada os ciclos já aí apre- 
sentados, referem-se alguns aspec- 
tos tecnológicos das instalações e é 
feita uma breve referência ao fun- 
cionamento de instalações a carga 
parcial. 

Um sector industrial muito ligado à 
termodinâmica aplicada é o ar con- 
dicionado. Assim, o Capítulo 2 
apresenta a teoria das misturas 
ideais e aplica essa teoria às mistu- 
ras de ar e vapor de água - o ar 
húmido. É definida a terminologia 
apropriada e apresentado o dia- 
grama psicrométrico. Numa se- 
gunda parte, são dados diversos 
processos com ar húmido, e os 
equipamentos onde ocorrem, aos 
quais são feitos balanços mássicos 
e energéticos. O condicionamento 
de ar é a aplicação final deste capi- 
tulo. 

A combustão é um fenómeno im- 
portante que a maioria dos enge- 
nheiros mecânicos virá a encontrar. 
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Assim, para a totalidade dos alunos 
o Capítulo 3 constitui uma introdu- 
ção à combustão. Começa por for- 
necer alguma informação sobre 
combustíveis e sobre a respectiva 
tecnologia de queima. Relembram- 
se algumas noções sobre reacções 
químicas e apresenta-se a análise 
energética da combustão, condu- 
zindo ao cálculo da temperatura de 
combustão adiabática. Finalmente 
considera-se o equilibrio quimico, 
as reacções de dissociação, e 
chega-se ao cálculo da temperatura 
de combustão adiabática conside- 
rando a dissociação dos produtos 
da combustão. 

A análise exergética (na terminolo- 
gia anglo-saxónica exergy analysis, 
second law analysis, ou availability 
analysis) é um método de análise 
que tem vindo a aumentar o seu 
campo de aplicação. As suas aplica- 
ções vao da simples análise de sis- 
temas energéticos, passando pelo 
estudo termo-económico de siste- 
mas e processos até ao planea- 
mento energético regional. O Capí- 
tulo 4 introduz o estudo dos 
sistemas energéticos em termos de 
exergia. Caracteriza a transferência 
de exergia associada ao trabalho e 
ao calor e a exergia associada a 
um escoamento. Define as compo- 
nentes cinética, potencial, física e 
química da exergia de um escoa- 
mento. Apresenta os processos irre- 
versíveis como correspondendo a 
uma perda (destruição) de exergia. 
Formula a equação de balanço de 
exergia e introduz o rendimento 
exergético de um sistema. Termina 
com alguns exemplos de aplicação 
da análise exergética a diversos 
processos e sistemas. 
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ABSTRACT 


Corrosion and scaling in geo- 
thermal systems are two of the 
more important problems that 
require the close attention of 
chemists, geologists, and ma- 
terials scientists. Of the many 
constituents in geothermal flu- 
ids, key chemical species are 
identified that account for most 
corrosion phenomena in geo- 
thermal systems. Based on 
analysis of actual geothermal 
materials test data and test 
methods, the performance of 
metals in geothermal fluids is 
presented for carbon and 
stainless steels, litanitum, nic- 
kel, copper, aluminium, and 
other alloys. 


RESUMO 


Os fenómenos de corrosão e 
de aparecimento de incrusta- 
ções nos sistemas geotérmi- 
cos constituem matéria de 
interesse para químicos, geó- 
logos e cientistas de materiais. 
Entre os constituintes dos flui- 
dos geotérmicos, identificam- 
-Se as espécies químicas mais 
susceptíveis de originarem 
problemas de corrosão nos 
meios geotêrmicos. Descreve- 
-se O comportamento de me- 
tais e ligas metálicas usadas 
nesses sistemas. 
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GEOTHERMAL FLUIDS 


The two geothermal fluids available 
are steam and water. Steam-domi- 
nated resources occur less fre- 
quentiy than water-dominated 
systems. Because of the greater 
abundance of water-dominated 
systems these resources will be 
emphasized in this paper. Hot dry 
rock and geopressured systems are 
mentioned only in passing. 


The water in a geothermal resource 
is of several types, as noted by Ellis 
and Mahon (1). Meteoric and mag- 
matic waters provide the primary 
sources and seem to be the most 
abundant. Meteoric water results 
from rain and snowball and differs 
from standard normal ocean water 
in its content of deuterium and 
heavy oxygen (O). Meteoric water 
generally has lower "O content 
than waters originating from mag- 
matic sources or meteoric water 
that has been in contact with silice- 
ous materials for long periods. The 
evidence to date suggests that geo- 
thermal fluids contain both meteoric 
and magmatic water; the relative 
amounts of each are characteristic 
of a particular resource. 


Consider now a drop of water falling 
through the atmosphere. The dissol- 
ved gases in this unit of meteoric 
water will be in equilibrium with the 
atmosphere as the water strikes the 
earth. The important chemical 
agents in this fluid are water, oxy- 
gen, and carbon dioxide; all are po- 
werful weathering agents. Water is 
both an active mineral solvent and 
a carrier for the dissolved minerals. 


Oxygen acts on sulfides and other 
minerals to form soluble products; it 
also acts on dead plant and animal 
matter to increase the rate of de- 
cay. Carbon dioxide is a mild oxidi- 
zing agent, but its primary function 
is found in its action on the kinds 
and amounts of carbonate species. 


Heturning briefly to plant and animal 
decay products, one finds that cer- 
tain of these organic substances 
are strong weathering agents. In 
particular, the humic and fulvic acid 
fractions are potent in briefing mine- 
ralis into solution and in forming 
relatively stable coordination com- 
pounds. Humic acids are characteri- 
zed by carboxvylic acid and phenolic 
functional groups. Singer and Na- 
vrot (2) found that a significant por- 
tion of the total metal content of 
metal-rich basalts was leached in a 
few hundred hours of contact with 
humic acid solutions at temperatu- 
res as low as 323 k. 


Within the near-surface environ- 
ment, the principal chemical wea- 
thering agents are water, oxygen, 
carbon dioxide, and decay products 
from plant and animal life. As these 
agents are transported deeper into 
the crust by a process termed elutri- 
ation, they bring about chemical 
changes and increase the load of 
dissolved solids and gases in the 
water and gases in the water gene- 
rally. 


As the water droplet penetrates still 
deeper into the crust, the tempera- 
ture of the surrounding roch increa- 
ses because of a heat flux from the 
molten core of the earth toward the 
cooler surface. This heat flux avera- 
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ges 0.063 J/m“ and results in a ther- 
mal gradient that varies with the 
thermal conductivity of the area. lts 
value is usually about 20 K/km. This 
heat flux results in heating of the 
meteoric water and in accelerating 
rates of reaction between water and 
rock. 


Magmatic water is released from 
solidifying magma. The heat source 
of the magma is the decay of radio- 
active elements, with uranium, tho- 
rium, and potassium providing most 
of the heat (3). As the radioactive 
elements are consumed, less decay 
heat is produced, and some of the 
magma crystallizes. There are local 
hot spots, as indicated by volcanic 
action, but for present purposes 
these can be ignored. A portion of 
the solidifying magma contains che- 
mically bound water. In water-domi- 
nated fields this magmatic water 
eventually cools enough to form a 
condensed phase. (In a steam-do- 
minated field a substantial portion of 
the geothermal fluid will remain as a 
vapor.) Although the liquid water will 
be free of oxygen and organic mat- 
ter, it will be saturated with dissol- 
ved rock matter and contain gases 
such as carbon dioxide and hydro- 
gen sulfide. The former is a product 
of carbonate decrepitation; the latter 
results from the decomposition of 
sulfide minerals. 


The hot meteoric water and the hot 
magmatic water will react with other 
rocks in the geological formation. 
As they pass from one stratum to 
another, some minerals will be de- 
posited and others dissolved. This 
geothermal fluid varies in composi- 
tion according to local lithology, so 
that geothermal fields show pro- 
found differences in water che- 
mistry. Furthermore, there are varia- 
tions from well to well within a field, 
and any single well may vary in wa- 
ter composition with time and with 
the rate of water production. 


Typically, chloride and sulfate are 
the most abundant anions, whereas 
sodium and calcium are the most 
abundant cations. Silica is the most 
abundant nonionic species in most 
geothermal waters. The ionic spe- 
cies are formed by the passing of 


the water through beds of evapori- 
tes and by the decomposition of va- 
rious minerals. The silica results 
from dissolution of a number of mi- 
neral species, including quartz. The 
ionic strength of most geothermal 
fluids is so high that serious deviati- 
ons from ideal models, such as the 
Debye-Huckel theory, have been 
observed. 


The hot “rock soup” is the fluid that 
must be worked with. This is the 
fluid from which heat is abstracted 
for electrical power production, for 
space heating, or for process heat 
in a number of applications. This 
paper addresses the character of 
this rock soup and materials for 
containing, transmitting, and han- 
dling it. 


SCALING AND CORROSION 


Once the geothermal fluid is availa- 
ble for use, two important technolo- 
gical problems confront us; these 
are scaling and corrosion. The sca- 
ling minerals are typically silica and 
calcite, although sulfides are impor- 
tant in some fields. Corrosion pro- 
ducts are a second source of scale. 
Both of these sources of scale are 
important because of the adverse 
effect of scale on heat transfer and 
pumping efficiencies. 


The species in the water that are of 
greatest interest in relation to corro- 
sion are hydrogen ion, chloride ion, 
hydrogen sulfide, carbon dioxide, 
oxygen, transition metal ions, am- 
monia, and sulfate ion (4). Some 
generalizations related to the corro- 
sive effects of these species are no- 
ted here: 


Hydrogen lon — The general corro- 
sion rate of carbon steels increases 
rapidly with increasing hydrogen ion 
(decreasing pH), especially below 
pH 7. The passivity of many alloys 
is pH dependent. Breakdown of 
passivity at local areas can lead to 
serious forms of attack, such as pit- 
ting, crevice corrosion, and stress 
corrosion cracking. 


Chloride — Chloride (Cl) causes lo- 


“cal breakdown of passive films, 


which protect many metals from 
uniform attack. Local penetration of 
this film can cause pitting, crevice 
corrosion, or stress corrosion crac- 
king. Uniform corrosion rates can 
also increase with increasing chlo- 
ride concentration, but this action is 
generally less serious than local 
forms of attack. 


Hydrogen Sultide —- Probably the 
most severe effect of hydrogen sul- 
fide (H,S) is its attack on certain 
copper and nickel alloys. These al- 
loys have performed well in seawa- 
ter but are practically unusable in 
geothermal fluids containing H,sS. 
The effect of H,S on iron-based ma- 
terials is less predictable. Accelera- 
ted attack occurs in some cases 
and inhibition in others. High- 
strength steels are often subject to 
sulfide stress cracking. Hydrogen 
sulfide may also cause hydrogen 
blistering of steels. Oxidation of H,S 
in aerated geothermal process stre- 
ams increases the acidity of the 
stream. 


Carbon Dioxide — In the acidic re- 
gion, carbon dioxide (CO,) can ac- 
celerate the uniform corrosion of 
carbon steels. The pH of geother- 
mal fluids and process streams is 
largely controlled by CO,. Carbona- 
tes and bicarbonates can display 
mild inhibitive effects. 


Ammonia —- Ammonia (NH,) can 
cause stress corrosion cracking of 
some copper alloys. It may also 
accelerate the uniform corrosion of 
mild steels. 


Sulfate — Sulfate (SO) plays a mi- 
nor role in most geothermal fluids. 
In some low-chloride streams, sul- 
fate will be the main aggressive 
anion. Even in this case, it rarely 
causes the same severe localized 
attack as chloride. 


Oxygen — The addition of part-per- 
billion quantities of oxygen (O.) to a 
high-temperature geothermal 
system can greatly increase the 
chance of severe localized corro- 
sion of normally resistant metais. 
The corrosion of carbon steels is 
sensitive to trace amounts of oxy- 
gen. 


Transition Metal lons — “Heavy” or 
transition metal ions might also be 
included as key species. Their ac- 
tion at low concentration on most 
construction materials is ill defined. 
However, the poor performance of 
aluminium alloys in geothermal flu- 
ids may be due in part to low levels 
of copper or mercury in these fluids. 
Some oxidized forms of transition 
metal ions (Fe*, Cu*?, and others) 
are corrosive, but these ions are 
present in the lowest oxidation state 
(most-reduced form) in geothermal 
fluids. 


Oxygen can convert Fe? to Fe*, 
which is another reason to exclude 
oxygen from geothermal streams. 


Corrosive attack of metals in geo- 
thermal systems can occur in seve- 
ral of the following forms (4 - 6): 


Uniform Corrosion — This is a gene- 
ral, allover attack of the metal sur- 
face. Uniform corrosion is often pro- 
moted by chloride, carbon dioxide, 
oxygen, or ammonia. 


Pitting — Pitting is a localized form 
of attack, which results in the deve- 
lopment of small pits in the metal 
surface. Pitting is often associated 
with the breakdown of a protective 
film or surface scale. 


Crevice Corrosion — Crevice corro- 
sion is similar to pitting in that itis a 
localized attack. Unlike most other 
forms of corrosion, it is geometry 
dependent and forms in the crevi- 
ces of equipment. 


Stree Corrosion Cracking — Stress 
corrosion cracking (SCC) is a catas- 
trophic type of failure promoted by a 
combination of tensile stress and 
the presence of chloride ion in the 
environment. The presence of oxy- 
gen and increasing temperature in- 
crease the severity of attack. 


Sulfide Stress Cracking — Sulfide 
stress cracking (SCC) is a catastro- 
phic failure that results from expo- 
sure under stress of susceptible 
materials to environments contai- 
ning H,S in an aqueous phase. In 
contrast to stress corrosion crac- 
king, SSC decreases in severity 
with increasing temperature, but 
oxygen may have little, if any, effect 
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on the SSC mechanism. On the 
other hand, low pH greatly accele- 
rates failure. 


Hydrogen Blistering — Rupture of 
metallic materials results when 
hydrogen trapped in voids accumu- 
lates at a sufficient pressure. The 
material need not be under stress 
for hydrogen blistering to oceur. 


Intergranular Corrosion — Intergra- 
nular corrosion is preferential corro- 
sion at or adjacent to grain bounda- 
ries, with little or no attack on the 
bodies of the grains. The alloy di- 
sintegrates (grain fall out) or loses 
its strength, or both. 


Galvanic Corrosion — This occurs 
when two metals are electrically 
connected. Corrosion of the less 
noble material will be accelerated. 
Materials may be ordered in a gal- 
vanic series (Dy increasing nobility) 
to help in materials selection. Care 
must be taken because the order of 
metals may change with variations 
in chemistry and temperature. 
Corrosion Fatigue — Corrosion fati- 
gue is premature fracture when 
cyclic stresses are imposed on a 
material in a corrosive environment. 
The corrosion fatigue limit is the 
greatest unit stress that may be ap- 
plied under given conditions of 
stress, rate of stress application, 
temperature, and corrosive environ- 
ment without causing the material to 
fail in a given number of cycles of 
stress (7). The combined effects of 
cyclic stress and corrosion are often 
far more severe than the simple 
sum of their actions. 


Exfoliation — This involves the for- 
mation of discrete layers of corro- 
sion products or of metal separated 
from the lattice by corrosion pro- 
ducts. The layers may break loose, 
damaging downstream componenis. 


MATERIALS SELECTION 


The selection of materials for appli- 
cations in geothermal environments 
is one area which chemists, metal- 
lurgists, and materials scientists can 
make a significant contribution to 
the geothermal effort. The respon- 
ses of plain carbon and low-alloy 


steels, stainless steels, nickel-based 
alloys, copper-base alloys, alumi- 
nium, titanium, and other materials 
are briefly reviewed in relation to 
common unflashed geothermal flu- 
ids. 

Plain Carbon and Low-Alloy Steels 
— Experience with these materials 
indicates that severe pitting and 
crevice corrosion occur, in addition 
to less severe uniform corrosion. 
The alloyed materials show more 
resistance but should not be selec- 
ted without a site-specific test that 
simulates the intended usage of the 
materials. These materials are sub- 
ject to hydrogen blistering, especi- 
ally in sulfide environments, and to 
sulfide stress cracking in the higher 
strength grades. Galvanic corrosion 
is strongly dependent on the cou- 
pling material, and, in our experi- 
ence, the corrosion of low-carbon 
steels is greater when the cathodic 
material is brass rather than stain- 
less steel. 


Stainless Steels — Ferritic alloys are 
subject to pitting and crevice corro- 
sion, especially in stagnant conditi- 
ons. However, some high-alloy ferri- 
tic stainless steels have shown very 
good resistance to a moderate-tem- 
perature geothermal fluid (8). Auste- 
nític alloys are prone to stress cor- 
rosion cracking, and at least one 
instance of this form of attack has 
been seen in geothermal usage. 
These alloys are also subject to lo- 
calized and intergranular corrosion. 
Testing of these alloys should in- 
clude the evaluation of oxygen in- 
trusion into the geothermal fluid. 
Martensitic alloys suffer from the 
same forms of corrosion as ferritic 
and austenitic stainless steels. 
Stainless steels should be tested in 
the proposed environment prior to 
being selected for use. Further, the 
testing should simulate the effects 
of oxygen on the alloys to evaluate 
their response to a stress corrosion 
cracking environment. 

Nickel Alloys — High-nickel alloys 
are frequently used to combat se- 
vere corrosion problems. The nic- 
kel-chromium-mobybdenum alloys 
appear to be the most applicable to 
high-temperature geothermal fluids. 
Inconel 625 and Hastelloy C-276 


are particularly resistant to corro- 
sion. These alloys can tolerate very 
high flow rates and occasional aera- 
tion. Similar alloys containing iron in 
place of molybdenum face competi- 
tion from the most resistant stain- 
less steels, but may find applicati- 
ons where their mechanical 
properties are desirable. Cupronic- 
kels will have limited usefulness in 
geothermal streams containing even 
trace quantities of H,S. 


Copper Alloys — The use of copper 
alloys in geothermal fluids is seve- 
rely limited by the relatively high 
concentrations of sulfide found in 
most sources. The Raft River, 
Idaho, KGRA, with a low sulfide 
concentration of 0.1 ppm, appears 
to be an exceptional case. Even in 
this fluid, the performance of cop- 
pernickel alloys has been very poor. 
Although some nickel-free brasses 
and bronzes have performed well in 
testing, dealloying of some copper 
alloys has been observed (8). 


The cracking of some copper-based 
alloys while exposed to ammonia or 
ammonia derivates is known as “se- 
ason cracking” (9) or ammonia 
SCC. Provided NH, and NH,* levels 
are low, cracking by ammonia spe- 
cies in the liquid phase of geother- 
mal fluid should be limited to stag- 
nant areas. Susceptibility of 
copper-zinc brasses to cracking in- 
creases with increasing zinc con- 
centration. 

Aluminium Alloys — Aluminium atl- 
loys have been recommended for 
desalination applications. Their use 
in geothermal environments cannot 
be recommended because of the 
very high corrosion rates that have 
been encountered. 

Titanium Alloys — Titanium and its 
alloys are among the most resistant 
materials tested to date. However, 
experience in other environments 
suggests that at temperatures 
above 573 K these materials are 
subject to crevice corrosion. Fluori- 
des specifically attack titanium, and, 
therefore, this material should be 
avoided where fluoride-rich waters 
are encountered. 


Other Alloys — Cobalt-base alloys 
are noted for their hardness and 


can be used in geothermal environ- 
ments when this characteristic is 
desirable. Zirconium has very good 
corrosion resistance in a number of 
environments, but the cost of these 
materials is such that no applicati- 
ons to geothermal environments are 
projected. 


CONCLUSION 


Corrosion and scaling in geothermal 
systems are two of the more impor- 
tant problems that require the close 
attention of chemists, geologists, 
and materials scientists. These pro- 
blems result from the nature of the 
geothermal fluids encountered as 
systems are designed and built to 
take advantage from these resour- 
ces. The geothermal fluid, in turn, is 
a result of the environment from 
which it is extracted. The history of 
the fluid as it comes into contact 
with various minerals is the key to 
understanding the tendency of 
these fluids to promote scaling and 
corrosion and, as a consequence, 
making rational materials selections 
for plant construction. Of the many 
constituents in geothermal fluids, 
seven key chemical species are 
identified that account for most cor- 
rosion phenomena in geothermal 
systems. These species are: oxy- 
gen, hydrogen, sulfide, carbon dio- 
xide, ammonia, chloride, sulfate and 
hydrogen ion concentrations. Based 
on analyses of actual geothermal 
materials test data and test me- 
thods, the performance of metals in 
geothermal media is presented for 
carbon and stainless steels, tita- 
nium, nickel, copper, and many 
other alloys. Eventually, sufficient 
data on and insight into the che- 
mistry of geothermal fluids will be 
available to permit adequate materi- 
als selection without longterm corro- 
sion tests. 
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ABSTRACT 


From the Takeno's study of the 
plane gravitational waves 
in general relativity, it Is 
observed that the cosmological 
term Y plays no role and 
hence its introduction in the 
field equations is not neces- 
sary. 


NIL CONTRIBUTION OF THE 
COSMOLOGICAL CONSTANT 
TO THE PLANE GRAVITATIONAL 
WAVES IN GENERAL 

RELATIVITY 


K. S. Adhav* and T. M. Karade 


Department of Mathematics, Nagpur University, Nagpur — 440010 — India 


INTRODUCTION 


The relativistic field equations first proposed by Einstein are 


where A is a cosmological constant (see Tolman 1934). 


In empty region with all the components of T, equal to zero we find that the 
above field equations reduce to the simpler equations 

R,=A&; (1) 
To the dimension of the solar system the effect of A term is negligible 
because the motions of the planets correspond with great precision to the 
field equations 

R,=0 (2) 

However, for regions of greater size a cosmological term may not neces- 
sarily be equal to zero. Therefore, it is worth to study the plane gravitatio- 
nal waves in reference to equations (1) instead of equations (2) with res- 
pect which Takeno (1961) has developed the theory of plane gravitational 
waves in general relativity. As a matter of fact, Takeno has shown that an 
empty space-time field equations (1) has no exact plane wave solution in 
the sense of the definition. We observe that this finding can also be dedu- 
ced right from violation of the signature convention of the underlying 
space-time. In this short communication we prove that the consideration of 
the plane gravitational waves in the empty region with equations (1) yields 
equations (2), thereby A contributing nil to the plane waves in general 
relativity. 


* Teacher fellow from the department of Mathematics, Shri Shivaji Science College, 
AMRBAVATI. 


33 


DINIL 


TEC) 


34 


FIELD EQUATIONS 


For generalisation we extend the Takeno's definition of a plane wave as follows: 
A plane wave g, is a non flat solution of equations (1) with the property that 


in some suitable coordinate system 


8; = Bi (2) , (X' =X, yzt) (3) 
which satisfy 
nú. Q 
g £,£,=D [i= 
' o x, (4) 
E = PA) s (Es PU SÃ£o O), (5) 
The signature convention is that of Takeno i.e. 
nú < O, (6) 
Ea tah > 0. (7) 
Oba Ebb 
Si &o É 
2) Bo Ba|<O 
e 4 23 Ea 
gu>0 , g=det(g)<o0 
(not summed fora and b; ab = 1,2,3). 
In the notations of Takeno, we get 
R " Es PP 
2M* 
(ab = 1,2) 
wir Pa=Esw wW=0g'+g": w=t+z6 
| N N 
A AD Liz—: Msg: M,=— a 
“-M “-M “M “-M 


M=wW-z) | N=W-zó 


and the bar (-) over a letter means the derivative with respect to Z. 


Then the field equations (1) for a,b = 1,2 become 


Pal Em 2 AN ME 


(8) 
which gives 
(p, ) =2nAM', (9) 
(p» ) =2A Mg» (10) 
and 
Ppp,=2nM'g, (11) 
1.8 PP e A M'(g,) 
Using equations (9) and (10) we get 
(2 A Mºg,) (2 A M'g.,) = a M(g, ) 
— &n&n= (815) 
Sim Bin] 0 
E Bo o» (12) 


which violates equation (7). 


Thus the study of the plane gravitational waves in the empty space-time 
field equations (1) ultimately reduces to equations (2) giving A = 0. 
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CONCLUSION 


It seems that the introduction of the cosmological term A in the field equati- 
ons of general relativity is not necessary as far as its role in the study of 
plane gravitational waves is concerned. 
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ABSTRACT 


In this paper, a eylindrically 
symmetric cosmological model 
is studied in the context of 
bimetric relativity, taking the 
source as a perfect fluid distri- 
bution. It is shown that the 
model does not exist in bime- 
tric relativity. 


A CILINDRICALLY 
SYMMETRIC COSMOLOGICAL 
MODEL IN BIMETRIC 
RELATIVITY 


Shailendra Deo 
Department of Mathematics, Nagpur University Campus, Nagpur-440 010, India. 


1 INTRODUCTION 


Rosen N [2] proposed a bimetric theory of relativity, where there exists two 
metric tensors at each point of space-time — g, which describes gravitation 
and background metric 7, which enters into the field equations directly with 
matter. Accordingly, at each space-time point one has two line elements 


ds* = g; dx' dx 
and do* = y, dx dx 


This theory is based on a simple form of Lagrangian and has a simpler 
mathematical structure than that of the general relativity. 


Here a cylindrically symmetric non-static Einstein-Rosen cosmological model 
is studied when the source is perfect fluid distribution. 


2 FIELD EQUATIONS AND MODEL 


Field equations of bimetric relativity formulated by Rosen N. [2] are 


| 
Ki=Ni > Nei =—8mkT; (1) 
where 

= | Cf gh 1 

N, E 2 F (g Bia da (2) 
A TR 

N=Nh. k=[E) (8) 

A 


a vertical bar (|) denotes the covariant differentiation with respect to y,. The 
energy-momentum tensor *, for a perfect fluid distribution is given by 


T;=(p+p) vv! + pg; (4) 


together with the orthogonality condition v,v* = -1, where v' is the four-velo- 
city vector of the fluid and p, ? are the proper pressure and the energy den- 
sity of matter respectively. 
Let us consider the cylindrically symmetric non-static Einstein-Rosen cosmo- 
logical model [1] given by 


d=e20P (JT? — dR?) — Re do? — eb do (5) 
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where «q. are functions of R and T, and the convention is 
x =Rx=4 W$=Zand x!=T. 
The flat-metric corresponding equation (5) is 
do” =dT* —dRº — R'dg” — dz” 


For this metric non-vanishing y-Christoffel symbols are 


2 3 


Using equations (1) to (7) the field equations yield 
Rº Sinh2a = — 8rkp 
o" +Ria- & -Rº Sinh2a = — 8rkp 


q +R'o- 4 -2pr-2RIB'+2 P =-Brkp 


O = Brkp 
where cíbi da di da 
dR OR - 
and 
o da o 
U=— , Q=- etc. 
dT JT 
Equation (11) implies 
p=U 


It shows that there is no contribution from perfect fluid distribution to the 
cylindrically symmetric non-static Einstein-Rosen cosmological model in 


bimetric relativity. The field equations (8) to (10) become - 
R“ Sinh2a = O 
o +R'o' - & —-Rº Sinh2a = 0 
e +R'o- O -25"-2RºB'+2 p = 0 


From equation (13) we get 


E 
| 
Õ 


Equation (15) and (16) imply 
Br+R'B'- É =0 
For solving equation (17) 


Case 1: Let us consider  B=Hírty)+G(t) 
then the equation (17) gives the Bessel equation 


H" + ty -H=0 
R 


representing cylindrical waves, provided 
6 =0 
The Bessel equation (18) results in 


H = a JAhR)Cos(ht + cj) + bY;(hRjSindht + 5) 


(6) 


(7) 


(8) 
(9) 


(10) 
(11) 


(12) 


(13) 
(14) 
(15) 


where Jo(hR) and Y,(hR) are the Bessel functions of first and 
second kind of order zero respectively, h is the frequency and a, b, 
C,, Ca are arbitrary constants. 


On integration, (19) gives 
G=ct+c (21) 
where c,, c, are arbitrary constants of integration 
Then 
B=a Jy(hR)Costht + e)+b Yo (hR)Sin(ht+c;)+est+e, (22) 


Case Il: Let us consider 5 = H(r.t) + g(r). 
In this case again we obtain equation (18) provided, 
g” + R'g' =0 (23) 
it gives us 
E = Es (24) 
where c, is an arbitrary constant of integration. 
Now the solution turns out to be 


B = a J,(hRjCos(ht + c)+b Y,(hR)Sin(ht+cs) + es (25) 


where constants could be found using initial as well as boundary 
conditions from the above equations. 
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ABSTRACT 


In the present paper a bro- 
adband microwave reflecto- 
metric system presently 
beeing developed for the 
tokamak ISTTOK is pre- 
sented. The diagnostic will 
measure the plasma den- 
sity profiles and fluctua- 
tions. The basic principles 
of density profile evaluation 
with microwave reflecto- 
metry are reviewed and the 
application of the theoreti- 
cal model to ISTTOK is 
discussed. Also, the speci- 
fic features of the system 
and its main characteristics 
and diagnostic capabilities 
are presented, 


A REFLECTOMETRIC SYSTEM 
FOR THE TOKAMAK ISTTOK' 


P. Varela, A. Silva, S. Vergamota, F. Serra 
M. E. Manso, J. A. C. Cabral* and C. A. F. Varandas 


Centro de Electrodinâmica, Instituto Supenor Técnico, 1096 Lisboa Codex 


1 INTRODUCTION 


The production of a plasma in a tokamak during a few seconds aims to pro- 
duce nuclear energy from thermonuclear controlled fusion. The rate of 
fusion reactions is fundamentally determined by the ion temperature T, and 
the Lawson parameter n,t. where n, is the density of reacting particles and 
Tt. is the ratio of the stored energy to the input power or energy confinement 
time. 

As the plasma is electrically neutral (n.= n) the knowledge of both n, and 
the density of impurities (given by spectroscopy) leads to the determination 
of the density of reacting particles. 


On the other hand, t. depends strongly on the spatial evolution of the elec- 
tronic density (peaked profiles lead to better confinement times than flat 
ones). In addition the energy confinement time depends on localized fluctua- 
tions of the plasma parameters related to various instabilities. An increase in 
the level of fluctuations can be correlated with a decrease in T.. 


The measurement of the electron density profile nr) as well as of the cha- 
racteristics of density fluctuations are therefore of fundamental importance to 
understand the physical processes involved in the fusion plasma. 


The two standard techniques which are currently used to probe the electro- 
nic density in high temperature plasmas are Transmission Microwave Inter- 
ferometry and Thomson Scattering; they lack however in general information 
about the edge region. Although this could be in part provided by Langmuir 
probes, a radial plasma region still remains where no data can be obtained 
by any of these diagnostics. In addition Thomson Scattering and Interfero- 
metry require multiple viewing chords to obtain spatial information about n,. 
The increasing restrictions in access to the machines limits the number of 
viewing chords that can be installed resulting in poor spatial resolution of the 
measurements. 


A technique which can meet the requirements of reduced access and good 
spatial resolution is Microwave Reflectometry. Basically, reflectometry is 
based on the propagation and reflection of an incident microwave beam in 
the inhomogeneous plasma. The frequency of the incident wave f determi- 
nes the value of the electronic density at the reflection layer. From the 
phase shift q(!) that waves with frequencies f between O and F undergoes in 
the plasma the position of each reflection layer can be inferred. 


* The “ISTTOK-Project” is supported by “Instituto Superior Técnico”, “Junta Nacional de 
Investigação Científica e Tecnológica” and the “European Atomic Energy Community 
(EURATOM)”. 


* Also belonging to “Centro de Electrodinâmica — Instituto Nacional de Investigação Cien- 
tífica”. 
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Here we describe the development of a reflectometric system for the toka- 
mak ISTTOK [1] with the possibility of broadband microwave sweeping and 
fixed frequency operation. The main purposes of the diagnostic are: (i) the 
measurement of density profiles; (ii) the study of density fluctuations in the 
gradient and edge regions of the plasma; (iii) experiments of radial correla- 
tion of the fluctuations parameters. 


2 THEORY 


2.1 Basic principles of microwave reflectometry 


Microwave Reflectometry is based on the total reflection experienced by a 
electromagnetic wave on the cutoff layer, where the refractive index vanis- 
hes. Propagation is perpendicular to the magnetic field lines (k | 6). The 
waves can be launched with their electric field either parallel (E Il B) or per- 
pendicular (E 1 B) to the magnetic field corresponding respectively to the 
Ordinary mode (O) and to the Extraordinary mode (X). The refractive index 
for these two modes, within a cold plasma approximation (v,. << c) and 
neglecting the motion of the ions (m, <<m,) and the damping effect of collisi- 
ons (f >> v,), can be written as; 


Ni=p- te (1) 
tr 
. n t- ; | 
West 2d 
E da À (2) 
“Sa 


Here f. and f., are the electron plasma frequency and the electron cyclotron 
frequency: 
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where e and m. are the charge and mass of the electron, €£, is the free 
space permittivity and B is the magnetic field.The corresponding cutoff fre- 
quencies are 


fo=f ay 
EE TAP | 
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where f, and f, are the low and high cutoff frequencies for the X mode. In 
the density range of interest, the lower cutoff corresponds to very low fre- 
quencies and therefore only the upper cutoff is considered. Two possibilities 


exist thus for wave reflection, corresponding to the O mode and upper X 
mode cutoffs. For the ISTTOK reflectometer we have selected the O mode 
wave propagation. One of the reasons for this choice is that N, depends 
only on n,, and the density profile evaluation is more straightforward. On the 
other hand, N, depends both on n. and B. This dependence on Bis a seri- 
ous drawback of the X mode since it makes reflectometry information 
dependable on another diagnostic (to measure B). Also, O mode reflection 
layers occur at lower frequencies compared with X mode which bring additi- 
onal practical advantages, since microwave components are less expensive 
at lower frequencies. In the following we will refer only to the O mode. 


Figure 1: Four different regions can be considered in the plasma. In region 1 f >> f,, and 
the wave travels at the speed of light. In region 2 f > f,, and the wave propaga- 
tes with the plasma refractive index Nfr). Region 3 is the reflecting region, 
where the refractive index decay to zero. In region 4 f< f,, and the wave is 
evanescent. 


The position of the reflecting layer r.(f) is basically determined from the 
phase shifts q(/) which the probing waves suffered in its propagation through 
the plasma. The phase shift for each frequency f can be found by matching 
an approximate WKB solution of the wave equation 


COLO Mi NtEnj=O (8) 
dr ce | 


outside the reflecting layer to the exact solution near the reflecting layer. Let 
us consider four different plasma regions as shown in figure 1. In region 
1 f>> f. and the electromagnetic wave travels with N(r) = 1 and velocity c. 
In this region the electric field of the wave takes the simple form 


E(r)= à,e"" + Beto (9) 


where k, is the wave propagation vector in the free space, k=27 f/c. In 
region 2 we have f > f.. The refractive index is always positive, N(r) < 1, and 
the wave can propagate through the plasma. The WKB approximation may 
be used to solve equation 8 if the inequality 


3 a. | Er <<1 (10) 
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can be applied. In that case the solution of equation 8 is given by 
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In region 3 the refractive index decay very quickly to zero and the WKB 
approximation is not valid there. Thus, an exact solution of equation 8 is 
needed. Defining V(r.) as a small region around the reflection point r., where 
the refractive index can be considered linear and making the variable 
change [2] 
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equation 8 can be written in the form 


E (ay o 
CEM E (y)=0 (13) 
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The solution of the above equation is 
E 00) = 0,400) + BB) (14) 


where A(x) and B(x) are the Airy functions. In region 4 the refractive index 
is always negative and the wave is evanescent. However in this region the 
WKB solution can be used and the solution of equation 8 is 
kf) | Nf dadr ho Mr dedo 
Biro DD df de? (15) 
| V Ni r) 


The constants a and / are determined by the continuity equations between 
adjacent regions. The exact solution is valid in the reflecting region, Jr. £,1.+ 
el, while the WKB solution is valid in the propagating region. Matching the 
two solutions, it can be shown [3] that 


dry 
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where q(f) is the phase shift due only to the propagation path in the plasma 
— from a (where n,= 0) to the reflecting layer r.(f). The term x/2 is related to 
the non-metallic character of the reflection in the plasma. 


The width of the reflecting layer, (2), can be estimated as the width of the 
region where the exact solution and the WKB solution differ by less than 
- 1%. In terms of x and for a parabolic profile (2) is estimated by 2 < xy < 5. 
The spatial resolution of the measurements correspond to the error in the 
location of the layer where the electric field vanishes which is given by the 
difference between 7=0 (WKB solution) and y=-2 (exact solution). 


2.2 Density profile evaluation 


For the case of O mode wave propagation the density at the reflecting layer 
for each probing frequency fis given by 


E, FR A 
| (17) 
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The location of the reflecting layer can be determined by analytically inver- 
ting equation (16) which yields: 


gato + og (18) 
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For the case of X mode propagation the refractive index is given by equa- 
tion (2). Due to its greater complexity it is not possible to analytically invert 
equation (16) and therefore a numerical invertion has to be performed in 
order to obtain r. 


2.3 Density fluctuations measurements 


Density fluctuations affect the refractive index thus the characteristics of the 
propagating wave, namely its electric field. The importance of the fluctuati- 
ons depends both on the local refractive index of the plasma and on the 
relation between the spatial scale length of the fluctuations A=27 / |k|l and 
the probing wavelength 4: 


(i) Small-scale density fluctuations — When A <A a significant part of the 
incident energy is scattered, due to the non-linear interaction between 
the incident wave (mo, k,) and the Fourier component (92, k ) of the oscil- 
lating modes perturbing the plasma, with frequencies m,=0, + Q and 
along directions defined by k=kt k 


(ii) Large-scale density fluctuations — When A >> 4 the scattering effects 
are negligible and phase modulations (4%) occurs which are mainly due 
to fluctuations localized in the reflection region. 


Several methods have been developed to analyse the reflectometer informa- 
tion concerning the fluctuations in the low frequency range (f < 20 KHz) 
associated with large-scale MHD activity. A method based on the harmonic 
contents of the detected signal has been proposed [4] and the contributions 
to the phase perturbations of the localization of the fluctuations were studied 
in detail, both in the reflecting layer (46,) and in the propagation region (46,). 
For typical fluctuations 46, >> 4q,. Based on phase measurements [5], a 
method has been proposed where the localization of the fluctuations is esti- 
mated by comparison of measured Ad(f) with values computed assuming 
perturbed density profiles. The standard model for a MHD density perturba- 
tion is used. Another method [6] based on the spectral analysis of the detec- 
ted signals, aimed at the evaluation of the fluctuations level and localiza- 
tion. Later, some other methods using WKB approximation and based on 
the phase shift modulation in inhoÔomogeneous plasma density profiles have 
been developed [7], [8]. A more accurate modeling of the scattered field was 
recently obtained using the full wave equation, [9]. 

More recently the problem of wave scattering by the fluctuations close to a 
cutoff layer and its application to fluctuations measurements have been 
studied [10]. For probing waves propagating along the plasma gradient, the 
model shows that the reflectometry diagnostic is local and sensitive only to 
fluctuations with low wavenumbers, k< (w/c)sL,'. Turbulence, generally 
with higher wavenumbers, lead to weak localization except in steep gradient 
density profiles. However if propagation in the k, is allowed good localization 
of turbulence, with k,=Kk,.-k, can be obtained. 


3 REFLECTOMETRIC SYSTEM 


The specific features of the system are discussed in the next para- 
graphs where the theoretical model is analysed in its application to the 
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tokamak ISTTOK. The main characteristics of the system are also pre- 
sented, 


Toroidal field 


Poloidal field 


waill 


Plasma 


Major radius + Minor ' R 
radius 


Figure 2: Geometry of a tokamak plasma. The coordinate system used as well as the 
poloidal and toroidal magnetic fields are also shown. 


3.1 ISTTOK 


Figure 2 shows the typical geometry of a tokamak plasma. It has a toroidal 
configuration characterized by a major radius A, and a minor radius a. Taking 
into account this geometric configuration we can use a toroidal coordinate 
system and therefore the spatial coordinates of a tokamak plasma can be 
defined by a toroidal angle 6, a poloidal angle 6, and a minor radius r. 


The magnetic field in a tokamak has two components: (i) a toroidal magnetic 
field B, produced by external coils that varies inversely with R=R,+r and is 
nearly independent of 6; (ii) a poloidal magnetic field B, produced by the 
plasma current |, which depends on r. The main parameters of the tokamak 
ISTTOK are summarized in Table 1. In tokamaks B, >> B. 


Parameter Value 


Minor radius, a 8.5 cm 


Major radius, A, 46 em 


Machine aspect ratio Dipo 


Toroidal magnetic field, B,, 0.45 — 0.6 T 
Plasma current, l, up to 8 KA 
Maximum density, Nam 0.6 — 1x 10 mº 


Electron temperature, T, 150 — 250 eV 


Table 1: ISTTOK main parameters. 


